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ABSTRACT 

We present a submillimetre survey of seven high-redshift galaxy clusters (0.64 < z < 1.0) 
using the Submillimetre Common-User Bolometer Array (SCUBA) at 850 and 450 jam. The 
targets, of similar richness and redshift, are selected from the Red-sequence Cluster Survey 
(RCS). We use this sample to investigate the apparent excess of submillimetre source counts 
in the direction of cluster fields compared to blank fields, as seen in the literature. The sample 
consists of three galaxy clusters that exhibit multiple optical arcs due to strong gravitational 
lensing, and a control group of four clusters with no apparent strong lensing. A tentative ex- 
cess of 2.7 a is seen in the number density of submillimetre sources within the lensing cluster 
fields compared to that in the control group. Ancillary observations at radio, mid-infrared, 
optical, and X-ray wavelengths allow for the identification of counterparts to many of the 
submillimetre luminous galaxies (SMGs), and provide improved astrometry and redshift con- 
straints. Utilizing photometric redshifts, we conclude that at least three of the galaxies within 
the lensing fields have redshifts consistent with the clusters and implied infrared luminosities 
of^ IO^^Lq. The existence of submillimetre cluster members may therefore be boosting 
source counts in the lensing cluster fields, which might be an effect of the dynamical state 
of those clusters. However, we find that the removal of potential cluster members from the 
counts analysis does not entirely eliminate the difference between the cluster samples. We 
also investigate possible occurrences of lensing between background submillimetre sources 
and lower-redshift optical galaxies, though further observations are required to make any con- 
clusive claims. Although the excess counts between the two cluster samples have not been 
unambiguously accounted for, these results warrant caution for interpreting submillimetre 
source counts in cluster fields and point source contamination for Sunyaev-Zel'dovich sur- 
veys. 

Key words: submillimetre: galaxies - galaxies: clusters: general - galaxies: high-redshift - 
galaxies: starburst - galaxies: evolution 
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1 INTRODUCTION 

1.1 Infrared-luminous Galaxies within Galaxy Cluster Fields 

The initial discovery of a population of dusty star-forming field 
galaxies with extreme luminosities in the high-redshift universe 



2 A.G. Noble et al 



with the Subminimetre Common-User Bolometer Array (SCUBA; 
e.g. |Smail et al.|1997|[Barger et al.|1998| [Hughes et al. 1998) has 
since spawned numerous submillimetre studies (e.g. Blain et^aTl 
T999b , Dunne et al. 2000 , Bales et al. 2000 , Best 2002 ; Cowie et al. 
2002, .C hapman et al.n2002nScott et al.,|2002nBorys et al..2003| 
Webb e t al. 2003a| |2005|TCoppin et al.|2006[|Knudsen et al.[2006t 
aimed at determining the cosmic star formation history and con- 
straining galaxy evolution models. Deep blank-field surveys with 
SCUBA (e.g. [Scott et al || 20021 [Borys et al.|[2003l [Webb et al.j 
[2003a[|Coppin et al.[20 06') of high-redshift submillimetre luminous 
galaxies (SMGs) have provided robust statistics to estimate source 
counts down to the 2 mJy level at 850 \xm. The lower end of the sub- 
millimetre flux density distribution has been probed by exploiting 
massive cluster lenses at low redshift (typically aiz < 0.4) to mag- 
nify the population of distant SMGs behind them (e.g.'Smail et ahj 
1997) [1998 2002; [C owie et al. |2002[ [Chapman et al. 2002, Knud^ 
sen et al.[2006| ), thereby reducing the source density on the sky and 
illuminating faint sources. These cluster surveys have led the way 
in understanding the properties of SMGs through numerous contri- 
butions. The first identification of radio counterparts ( jlvison et"ar| 
p^98l [2000l |Smail et al.[2000) secured the astrometry of individual 
SMGs, and thus facilitated the first spectroscopic redshifts ( [Ivisonj 
[et al.|1998t,2000„Barger et sl\.\1999) . Moreover, follow-up observa- 
tions of SMGs from cluster surveys yielded the first CO detections 
([Prayer et al. |1998| [T9991 jlvison et stl'lOO ^ and the discov ery of 
extremely red object (ERO) counterparts ( Small et al.|1999| ). 

An elucidating comparison can be made between source 
counts in the field and in cluster environs. Consistently, the clus- 
ter counts lie above the blank-field counts, which could be the 
result of unaccounted-for lensing or infrared galaxies associated 
with the clusters themselves ( Borys et al.|2003[|Webb et al.|2005] ). 
This holds despite the fact that lower redshift cluster surveys typi- 
cally remove central dominant (cD) galaxies from counts analyses 
( [Cowie et al.|2002[|Smail et al.| 2002), or expect few cluster mem- 
ber detections ([Chapman et al. 2002 , Knu dsen et al. 2006 ). Clus- 
ter counts in single pointings or small maps can carry a potential 
bias due to human intervention in two ways. For instance, the ten- 
dency to not publish blank maps could bias the available sample 
of clusters towards more crowded fields. Secondly, real-time ob- 
servers are more inclined to increase the integration time on areas 
that look more dense in small maps, and therefore, introduce addi- 
tional sources that would otherwise go undetected. However, it is 
impossible to quantify these possible biases in the literature (but see 
§4.3| for a discussion with our own work); we will therefore assume 
they are not issues in the above surveys and instead focus on the 
two aforementioned physical reasons for enhanced cluster counts. 

With regards to the unaccounted lensing scenario, higher red- 
shift clusters are generally less efficient lenses (i.e. have a lower 
lensing magnification) than their low-redshift counterparts ( |B artel- 1 
[mann et aL]|1998| ). This implies that unaccounted-for lensing is 
less likely to produce higher counts in high-redshift clusters fields. 
However, high-redshift cluster surveys that have calculated negligi- 
ble lensing corrections for the SMGs within the fields still observe a 
tentative excess of source counts compared to the field ( ^Best|2002[ 
[Webb et al.|2005 t this work), though this relies on small number 
statistics. 

If the excess counts are instead due to submillimetre bright 
cluster members, we would expect the effect to become more pro- 
nounced in higher redshift clusters (i.e. the surveys of [Best|2002] 
[Webb et al.|r2 005, this work) given the following trends with in- 
creasing redshift: a rise in the fraction of star-forming galaxies (e.g. 
[Butcher & Oemler[1978[[Saintonge et al.,2008„Haines et al.[2009j ); 



an increase in field infall (e.g. Bower|1991[rEllingson et al.|200T| 
[Loh etal.[2 00¥) and younger inf ailing galaxies; and the less relaxed 
state of younger clusters. Clarifying the source of excess counts 
will not only improve the submillimetre source counts at the faint 
end, but also provide information on the extreme phases of galaxy 
evolution in dense cluster environments, and aid cluster cosmology 
studies, which require constraints on the point source contamina- 
tion within cluster fields, specifically Sunyaev-Zel'dovich surveys. 



1.2 Strong Lensing in Clusters 

The efficacy of a cluster to produce multiple strong arcs primarily 
depends on the cluster redshift, the central surface mass density, the 
distribution of background sources, and the symmetry of the struc- 
ture; these requirements generally imply that high-redshift clusters 
are poor gravitational lenses compared to their low-redshift coun- 
terparts. The critical surface mass density required depends on red- 
shift, with low-redshift clusters being the most efficient lenses for 
background sources at z ^ 1 ( |Bartelmann et al.^1998) . The fourth 
efficiency factor dictates that asymmetric clusters with enhanced 
substructure typically produce larger arcs compared to those from 
spherically symmetric cluster models ([ Bartelmann et al.|1995j ). 

In contradiction to this redshift trend, a subset of high-redshift 
clusters with multiple giant arcs has been discovered with the first 
Red-sequence Cluster Survey (RCS-1; [Gladders & Yee][2005] ) at 
z > 0.64, while at lower redshifts, the survey failed to detect any 
lensing clusters, implying the property may be linked to the pro- 
cess of cluster formation |Gladders et al.|2003] ). While the excess 
of high redshift super-lenses can be explained with a higher value 
of (78 ( |Li et al.|2006j ), the same effect can be produced through a 
population of massive, concentrated clusters ( Dala i et al.|2004| ) or 
clusters undergoing a major merger, with substantial substructure. 
In each scenario, the central lensing cross-section can be temporar- 
ily enhanced, thereby boosting lensing efficiency (e.g. [Torri et aL] 
2004 ). Thus, these super-lenses may represent clusters which are 
dynamically unrelaxed or merging, both of which have been shown 
to increase the activity within their cluster galaxies ( iBlain et alT] 
[1999a||OweR et al.| 1999 [[Miller & Owenl2003l ). 

1.3 The Study 

To investigate the possible excess of SMGs in cluster environments 
and the presence of giant optical arcs from strong lensing effects, 
we analyse seven galaxy clusters of similar richness and redshift 
(0.640 < z < 1.045) drawn from the RCS-1 sample (see ^TTJ. The 
clusters are divided into two subsamples: the strong lensing sample 
and the control sample. The former consists of three galaxy clusters 
that exhibit a high incidence of optical arcs due to strong gravita- 
tional lensing within the central ^ 30 arcsec compared to theoret- 
ical predictions within a ACDM cosmology framework ( [Gladders] 
[et al. 2003 ). As a means of assessing the importance of the the 
cluster lensing efficiency, we also include four clusters (the control 
group) that cover the same range in richness and redshift, but lack 
any evidence of strong lensing effects. The original goal of the sur- 
vey was to obtain five clusters for each sample, but SCUBA ceased 
working before the remaining observations were completed. The 
entire sample is listed in Table [T] along with relevant cluster prop- 
erties. 

If excess source counts in cluster fields derive from lensing ef- 
fects due to the cluster potential, we should not observe this effect 
in our seven cluster fields, assuming our sample contains inherently 
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Table 1. A summary of the seven galaxy clusters of the submillimetre survey. 



Cluster Name 


Subsample 


Redshift" 


£)gc 

(Mpc^-'^'^) 


(1014 M0) 


(keV) 


~3 cr Depth^ 
(mJy/beam) 


Map Area-^ 
(arcmin^) 


Ancillary Data^ 


RCS 022434-0002.5 


Lensing 


0.773 


782 ± 201 


3.3 + 1.3 


^ 0+1-2 


3.1 


7.02 


SMlRXz'RcVB 


RCS 141910+5326.2 


Lensing 


0.640 


1400 ± 256 


9.5 + 3.3 


A r+0.4 

4-5Io.3 


4.3 


6.64 


SMlRXz'Rc 


RCS 231953+0038.0 


Lensing 


0.9024 


2107 ± 433 


18.7 + 8.3 


6 2+°-9 


3.1 


6.75 


SMlRXz'RcVB 


RCS 112225+2422.9 


Control 


0.940* 


1312 + 350 


7.1 + 3.5 




2.8 


6.65 


SMlz'RcVB 


RCS 132629+2903.1 


Control 


1.045* 


1901 +451 


13.8 + 6.9 


1 6+°-'^ 

-^•"-0.3 


2.5 


6.81 


SMlXz'Rc 


RCS 215248-0609.4 


Control 


0.650 


935 + 241 


4.7+1.9 




3.0 


6.62 


SMlz'RcVB 


RCS 231831+0034.2 


Control 


0.7518 


1153 + 252 


6.2 + 2.4 


6 1+1-^ 


2.7 


6.60 


SRXz'Rc 



"Redshifts are all spectroscopic except for those denoted with an asterisk (*), which are based on the location of the red sequence and accurate to 
10 per cent. 

^The amplitude of the cluster-centre galaxy correlation function. 

^The mass contained within the radius at which the cluster mass density is 200 times larger than the critical mass density of the universe. The value 
of M200 is calculated from th e velocity dispersion , which is correlated with Bgc ( |Yee & Ellingson|2003) . 
^The X-ray temperature, from |Hicks et al.||2008) , for the five clusters with X-ray coverage. 
^The average 3 a depth over the 850 |^m maps. 

^ The usable area of the 850 |xm maps over which the noise the uniform. 

^The ancillary data coverage for each field (see §2.3) . The symbols are as follows: SCUBA 850- and 450-|Lim imaging (S); MIPS 24-|Lim imaging 
(M); four-channel IRAC coverage (I); 1.4-GHz imaging with the VLA (R); Chandra X-ray data (X); original RCS filters (z' , Re); optical filters 
(V, B). Photometric redshifts are only available for clusters with all four optical filters. 



bad lenses compared to local clusters. Although seemingly contra- 
dictory to our claim above, the efficiency of our lensing sample to 
produce arcs arises within the central ^ 30 arcsec of the cluster, not 
over the entire field. The count excess in this study, however, is seen 
over a larger area (^ 7 arcmin^), where the strong lensing boosting 
factor is expected to be quite small. Certainly, we should not see an 
excess in the control fields, which show no signs of strong lensing. 

Alternatively, if a trend towards higher counts exists in all 
seven clusters, then it is likely an evolutionary effect with cluster 
galaxy populations. Moreover, if the excess is greater or limited to 
the super-lenses, then perhaps the physical mechanism driving this 
effect is the dynamical state of the cluster or excess lensing due to 
the cluster potential (but only in the cluster centre). We can there- 
fore attempt to disentangle lensing and cluster contamination, and 
further isolate the driver of contamination in the latter case. 

We describe our observations and data reduction for the 
submillimetre and the ancillary data sets in ^ In Q we out- 
line the submillimetre data analysis, and present 850 ^tm source 
counts in ^ The counterpart identification process with the multi- 
wavelength follow-up data is presented in ^ We discuss both 
the candidate submillimetre cluster members and potential cases 
of lensed SMGs, and their respective effects on the submillimetre 
source counts in ^ We conclude in ^ and include a discussion of 
each SMG in Appendix [A] We assume a flat, Qm = 0.3, Ha = 0.7 
cosmology, with Hq — 70kms~^ Mpc~^ throughout the paper. 



2 OBSERVATIONS AND DATA REDUCTION 

2.1 The Red-sequence Cluster Survey Sample 

To study the above issues, we have observed, using SCUBA at 850 
and 450 |am, a unique sample of seven galaxy clusters, drawn from 
the RCS-1 ( Gladders & Yee|2 005 ). The three lensing clusters (RCS 
0224-0002, RCS 1419+5326, and RCS 2319+0038) all lie within 



0.640 < z < 0.9024 and contain at least two optical arcs from 
strong-lensing. The remaining four clusters (RCS 1 122+2422, RCS 
1326+2903, RCS 2152-0609, and RCS 2318+0034) have redshifts 
within 0.650 < z < 1.045 and are used as a control sample. 

The RCS utilizes the well-defined relation between colour and 
magnitude of early- type galaxies, known as the red- sequence, to 
find clusters. The method is robust with a contamination rate of < 5 
per cent and provides a redshift estimate to ^ 10 per cent ( |Glad-| 
|ders & Yee(2000l|Gilbank et al.|2007) . Cluster significance relies 
on detecting simultaneous overdensities in colour, magnitude, and 
position. 

Launching the largest systematic search for galaxy clusters up 
to 2; ^ 1 at the time, the RCS-1 imaged ^ 92 deg^ and discovered 
^ 1000 galaxy clusters using only two optical filters {Rc and z'). 
The bands span the 4000 A break and provide superior colour dis- 
crimination at high redshift (up to 2; ^ 1) compared to the V and 
/ filter pair typically used in optical surveys. Moreover, the sur- 
vey includes an optical richness estimate based on the Bgc param- 
eter ( [Longair & Seldner|1979| ), the amplitude of the cluster-centre 
galaxy correlation function, which has been shown to be a good 
proxy for mass ( |Yee & L6pez-Cruz|1999||Yee & Ellingson|2003| ). 

2.2 Submillimetre Data 

2.2.1 SCUBA Observations 

The seven clusters were observed on the James Clerk Maxwell 
Telescope (JCMT) using SCUBA, at 850 and 450 |^m, over the 
course of several observing runs in 2001 and 2002. The data were 
reduced following the standard SCUBA User Reduction Facility 
(SURF) routine { [Jenness & Lightfoot|[T998| ). A 30-arcsec chop 
throw in right ascension was used for each cluster observation. 
Each map consists of either one or two overlapping pointings, re- 
sulting in areas all greater than ^6.5 arcmin^. The edges of each 
map, however, are unacceptably noisy due to lower integration 
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times and are not used in further analysis (see §2.2.2| ). Pointing sta- 
bility was maintained and checked after each major telescope slew. 
The sky opacity was monitored regularly through "skydips" every 
1-2 hours, allowing for a determination of signal attenuation due to 
a lack of atmospheric transmission. The opacity at 850 \xm ranged 
from roughly 0.15-0.29, with a typical value below 0.2. 

The presence of a 16-sample noise spike complicated the data 
reduction for the 2002 observations. This strongly correlated noise 
signal was removed from corrupted bolometers as described in 
[Webb et ar] ( |2005 1). The final corrected and uncorrected data were 
compared to the noncorrupted 2001 data and it was concluded 
that no robust point sources were either lost or introduced into the 
map, though the overall noise was reduced by ^ 20 per cent. The 
SHADES survey was similarly affected and also found the spike to 
have a negligible effect on the data ( [Coppin et al.|2006| ). 



2.2.2 850 iim Source Detection: Noise Analysis and Cleaning 

Robust source detection and minimization of spurious sources rely 
on accurate noise analysis. Due to the inconsistent noise proper- 
ties of SCUBA data and the generally low signal-to-noise (S/N) at 
which one is working, standard chop and nod techniques are not 
sufficient to remove all the atmospheric noise. Maps characterizing 
the noise are created from a set of Monte Carlo simulations for each 
bolometer time stream and a bootstrap algorithm (outlined in |Eales| 
let al.|2000| |Webb et al .|2003ll|2005] ). 

The jiggle-map scan pattern of SCUBA pointings produces 
uneven noise characteristics across the map; specifically, the edges 
of the map are much noisier compared to the centre due to the 
shorter integration time where the pointings do not overlap. To 
prevent falsely classifying spurious noise spikes as real sources, a 
mask is created to eliminate map regions from subsequent analysis 
where the noise exceeds 2mJy, primarily in the outer 15 arcsec of 
the map. Using maps with more uniform noise and even coverage 
helps maintain consistency throughout the data analysis process. 
The resulting map areas of the seven clusters range between 6.60 
and 7.02 arcmin^ (see Table [!}. It is within these areas that both 
source extraction and Monte Carlo simulations are performed. 

We create an initial source list of S/N > 3 detections from the 
uniform noise map area after convolving the raw signal map with 
a beam template, which consists of a 14-arcsec point source kernel 
and the negative off- source chops on either side. In low signal-to- 
noise images, a template-convolved map improves the accuracy of 
photometric parameters fPhillipp s & Davie s'1991), helps discrim- 
inate between real and spurious sources, and increases the signal- 
to-noise due to the inclusion of signal from the negative beams. 

We use an iterative cleaning process (see ,Eales et al.||2000[ 
|Webb et al.|2005l for more details) to isolate each potential source 
on the raw map, and separately convolve it with the beam template 
to reduce source confusion. Our final catalogue, shown in Table [2] 
consists of 26 sources above 3 a, with 2 — 6 SMGs per cluster field. 



2.2.3 450 iim Detection 

A unique advantage of SCUBA is its ability to simultaneously ob- 
serve two submillimetre wavelengths using a dichroic beam split- 
ter and two separate bolometer arrays, thus providing 450-|j,m 
data for free. In principle, the smaller beam associated with the 
shorter wavelength should result in improved positional accuracy of 
sources; strong atmospheric emission at 450 ^im, however, causes 
deterioration of data quality. Furthermore, telescope sensitivity to 



slight temperature changes cause high variability in the 450 ^im 
beam, thus complicating beam fitting. We therefore resort to simple 
smoothing techniques, rather than convolving the raw data with a 
beam template, when performing source detection. 

Despite this setback, our cluster sample contains fairly high 
quality data (i.e. less noisy than typical) at 450 \im, allowing source 
extraction and counterpart identification to the 850 ^tm emission. 
We adopt a search radius of 12 arcsec around an 850 ^im peak and a 
S/N cutoff of 2.5 cr at 450 \im to locate any 450 |^m peak emission 
associated with the 850 M,m SMGs. In the RCS 0224 field, Webb] 
|et al.| ( |2005| ) conclude that simulations of the 450 \xm positional 
uncertainties agree with those at 850 |^m (see §3.4| ), calculating a 
13 arcsec offset for the 95th percentile. The large search area of 144 
arcsec^, therefore, accounts for errors on positions at both wave- 
lengths. The low detection limit (2.5 a) ensures that we will not 
miss any 450 \im counterpart emission. 

For the 26 SMGs, we find 450 \im counterpart emission for 
12 sources, 10 of which exceed a S/N of 3.0 at 450 \im. Due to the 
smaller 450 ^im field of view (FOV), four SMGs fall outside the 
shorter wavelength coverage, and seven lie on the noisy edge of the 
field (see Table [21. We plot 450 M,m S/N contours on the 850 \im 
S/N map in Fig.[^(lensing clusters) and Fig. |2] (control clusters). 

We compute the probability of detecting random 450 jam 
emission within 12 ar csec of the SMGs using the P- statistic from 
)wnes et al. ( 1986| l, P — 1— exp(— Trnr^), where n is the sur- 
face density of 450 |j,m sources that have higher S/N than the can- 
didate identifications. In most of our cluster fields, we do not find 
any random detections with the required S/N. In RCS 0224, 1419, 
and 2318, we calculate P-statistics of 0.05, 0.09, and 0.04, respec- 
tively. This confirms that only 0-1 450 \xm candidate identifications 
per cluster are random associations. 



2.3 Ancillary Data 

In addition to the SCUBA imaging, we have ancillary observations 
across the spectrum, including optical, mid-IR, radio, and X-ray. 
Unfortunately, the coverage is not homogeneous over all the clus- 
ters. A summary of the multi- wavelength data available for each 
cluster is provided in Table [T] and the observations are described 
below. 



2.3.1 Optical Data 

Multi-band optical imaging for all clusters exists from the RCS sur- 
vey with the CFH12K camera on the Canada-France-Hawaii Tele- 
scope (CFHT) and the Mosiac II camera on the Cerro Tololo Inter- 
American Observatory (CTIO) Blanco telescope in z' and Rc fil- 
ters (Gladders & Yee 2005). Follow-up observations in V and B 



D5l 
-a(|F 



using the CFH12K camera ( [Hsieh et alT]2 005 ) have provided four- 
colour photometry for several optical sources within four cluster 
fields (RCS 0224, 1122, 2152, and 2319). The data were used to 
calculate photometric redshifts. For details on the observations and 
the photometric redshift fitting technique, see , Hsieh et al. ( 2005] ). 

2.3.2 Spitzer MIPS Imaging 

The Multiband Imaging Photometer for Spitzer (MIPS) at 24 |j,m 
has optimal capabilities to detect any mid-IR emission associ- 
ated with the 850 |j,m SCUBA sources. With a 5.4arcmin^ field 
of view and point spread function of ^ 6 arcsec, MIPS provides 
improved source positions compared to the SCUBA maps with a 
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Table 2. Positions and flux densities of the 850 \im submillimetre sources. 



R.A Dec. (measured) (deboosted) 850 \im (measured) 450 |xm 

Object Name" (J2000.0) (J2000.0) (mJy) (mJy) S/N (mJy) S/N 



SMM 0224.1 02 24 34.19 -00 03 25 6.2 6.2 6.1 14.4 3.7 

SMM 0224.2 02 24 33.46 -00 03 56 4.9 3.8 4.3 <54 

SMM 0224.3 02 24 28.06 -00 03 18 4.3 3.4 3.5 <115 

SMM 0224.4 02 24 29.79 -00 02 59 4.2 3.2 3.9 15.8 3.3 

SMM 0224.5 02 24 32.66 -00 03 42 3.3 2.4 3.4 18.7 2.8 

SMM 0224.6 02 24 32.93 -00 02 32 3.1 2.2 3.7 13.2 3.1 



SMM 1122.1 1122 29.42 +24 21 17 4.0 3.1 3.2 — — 

SMM 1122.2 11 22 31.69 +24 22 09 3.6 2.6 4.3 — — 

SMM 1122.3 1122 27.08 +24 22 41 2.5 1.5 3.4 16.0 5.6 



SMM 1326.1 13 26 28.18 +29 03 13 4.6 3.5 6.5 22.4 6.2 

SMM 1326.2 13 26 33.06 +29 02 20 3.5 2.5 3.6 <57 

SMM 1326.3 13 26 33.60 +29 03 23 2.9 1.7 3.4 <31 

SMM 1326.4 13 26 25.28 +29 02 41 2.5 1.5 3.2 <10 



SMM 1419.1 14 19 15.47 +53 24 55 8.3 8.3 5.2 — — 

SMM 1419.2 14 19 08.76 +53 24 56 7.3 7.3 3.0 — — 

SMM 1419.3 14 19 20.17 +53 26 29 5.0 4.5 4.5 19.3 2.7 

SMM 1419.4 14 19 05.51 +53 26 20 4.7 3.6 3.8 <16 

SMM 1419.5 14 19 14.69 +53 25 21 4.4 3.4 5.9 <13 



SMM 2152.1 21 52 46.62 -06 08 49 3.4 2.4 3.5 <26 

SMM 2152.2 2152 43.74 -0610 20 3.0 2.1 3.0 <54 



SMM 2318.1 23 18 27.89 +00 34 55 14.9 14.9 17.0 36.5 7.1 

SMM 2318.2 23 18 28.03 +00 34 18 5.0 4.5 5.6 14.5 3.0 

SMM 2318.3 23 18 32.63 +00 34 22 3.6 2.6 4.0 14.2 3.3 



SMM 2319.1 23 19 49.73 +00 37 57 7.6 7.6 5.5 26.1 4.4 

SMM 2319.2 23 19 56.86 +00 37 21 7.4 7.4 7.5 44.8 4.4 

SMM 2319.3 23 19 53.26 +00 38 09 3.3 2.4 4.7 <12 



"Sources in italics are detected below 3.5cr at 850 \xm and are therefore less reliable. In particular, we expect four of our sources to be spurious 
(see |33}. 

^The measured 850 |j.m flux, as described in ^2.2.2] 

^The deboosted 850 |Lim flux, which is a more accurate representation of the source (see §3.3) . 

^Objects without a 450 |j.m detection and those lying close to the edge of the usable 450 jiim field (where the noise is high) are given a 3 cr 
upper limit on the 450 |j.m flux. The four objects that lie completely off the 450 |^m field of view are denoted by a soHd dash ( — ). 



15 arcsec beam. These wavelengths typically trace the same popu- 
lations since both probe dusty star formation. 

MIPS 24 \im observations on the Spitzer Space Telescope 
(proposal ID 30940) have yielded complementary data for six of 
the seven SCUBA clusters (RCS 2318 was inadvertently left out 
of the observing list), with integration times ranging from 2000- 
3300 seconds per pixel. These data are discussed in more detail 
in Webb et al. (in preparation). Briefly, the MIPS images were re- 
duced using a combination of the Spitzer Science Center's MOPEX 
software and a suite of IDL routines we developed to further opti- 
mize background subtraction. Source catalogues of the entire MIPS 
field contain flux densities estimated with DAOPHOT, the Domin- 
ion Astrophysical Observatory's stellar photometry package ( Stet-| 
|son|198 7\ and positions using a photometry pipeline ( Y ee|1991| l. 
Monte Carlo simulations indicate 80 per cent completeness limits 
of ^ 70 pJy, or Lir 10^^ Lq, for all clusters. 



2.3.3 Spitzer IRAC Imaging 

The Infrared Array Camera (IRAC) on Spitzer operates at 3.6, 4.5, 
5.8, and 8.0 ^im. We obtained IRAC data (program ID 20754) in 



all four channels for each SCUBA cluster, except RCS 2318, with 
120 second integrations for all clusters. The images were processed 
using the Spitzer Science Center's IRACproc software suite ( |Schus-| 
Iter et al.||2006t . The source finding was completed with a simi- 
lar source detection and photometry pipeline that was used for the 
MIPS catalogues. The images reach a 3 a depth of ^ 4.2, 6.0, 27, 
26 |j,Jy in channels 1-4, respectively. 



2.3.4 Radio Imaging 

Deep 1.4 GHz pseudo-continuum (L-band) imaging from the Na- 
tional Radio Astronomy Observatory's (NRAOjjVery Large Array 
(VLA) of the three lensing clusters exists from observations in 2003 
through 2005. Using the 'A' and 'B' configurations we achieve a 
resolution of ^ 4.5 arcsec at 1.4 GHz and reach down to a noise 
level of ^ 15 ^tJybeam"^ after approximately 16 hours of integra- 
tion time on each cluster. 



^ NRAO is operated by the Associated Universities, Inc., under a coopera- 
tive agreement with the National Science Foundation. 
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Figure 1. 850 |Ltm S/N maps (convolved with the beam template) for the three strong lensing clusters. Overlaid are 450 |Lim S/N contours (smoothed with an 
8 arcsec Gaussian) in steps of 0.5, beginning with 2.5 a (the significance limit used for 450 \im counterpart emission), and steps of 1.0 above 4.0 a. Note the 
smaller field of view for 450 |Ltm, evident from the large orange borders. Significant 850 |j.m detections (within the usable map area — see §2.2.2| a re evident as 
light blue or white areas (i.e. ^ 3.0 cr) flanked by negative (black) lobes, and marked with numbers corresponding to their SMM IDs in Table |2] The location 
of the brightest cluster galaxy (BCG) is marked by an 'x' in each field; note there are multiple BCGs in some clusters. Clockwise from the upper left corner, 
the clusters are RCS 0224, 1419, and 2319. In each map, north points up and east is towards the left; the map centres are given by the coordinates in the cluster 
names of Table[T] 



Target data were recorded every 5 seconds and periodically in- 
terrupted with observations of a nearby phase target in order to ob- 
tain calibration for phase, amplitude and bandpass. Flux calibration 
was performed using standard unresolved sources with known flux 
densities; 3C48 was observed for RCS 0224 and 2319, while 3C286 
was used for RCS 1419. 

Calibration and cleaning/imaging were completed using stan- 
dard tasks in the Astronomical Imaging Processing System (AIPS) 
software. The resulting clean images are ^ 34 x 34 arcmin^ (for the 
central field), with RCS 1419 and RCS 2319 having slightly larger 
areas at each side due to mosiacing. The RCS 2319 image also 
extends into the RCS 2318 field; the RCS 2319 radio catalogue 
therefore includes sources within the two SCUBA fields. 

Source extraction and flux determination are achieved with the 
AIPS task 'SAD' (Search And Destroy). To reduce the detection of 
spurious sources, we first generate a S/N map from the clean image. 
With the original clean image and noise map as inputs, the 'SAD' 
task locates potential sources down to a specified level and uses a 
Gaussian fitting algorithm to determine the least square Gaussian fit 
to the source. In the case of extended sources, up to four Gaussian 



components can be fit and integrated flux densities are summed. 
We search down to a S/N of 3.5 cr for each cluster. 

2.3.5 X-ray Data 

X-ray observations, described in |Hicks et aT] p008| ), with Chan- 
dra's Advanced CCD Imaging Spectrometer (ACIS) were taken 
between 2002 and 2005 of five (RCS 0224, 1326, 1419, 2318, and 
2319) of the seven clusters with individual exposures ranging from 
50 — 100 ks. After initial cleaning, 0.3 — 0.7 keV flux images were 
created of the central 4x4 arcmin^ of each cluster. The point source 
sensitivities range from 2.5 — 5.1 x 10~^^ ergs s~^ cm~^. 

3 SUBMILLIMETRE DATA ANALYSIS AND RESULTS 
3.1 Monte Carlo Simulations 

Monte Carlo simulations allow for determination of the uncer- 
tainties on source properties, such as flux densities and positions 
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Figure 2. 850 |Ltm S/N maps (convolved with the beam template) for the four control clusters with the same explanation as Fig.[T] Clockwise from the upper 
left corner, the clusters are RCS 1 122, 1326, 2152, and 2318. Note that we have not been able to identify the BCG for RCS 2318. In each map, north points up 
and east is towards the left; the map centres are given by the coordinates in the cluster names of Table[T] 



(see §3.3| and §3.4[ ), and signal map characteristics, such as com- 
pleteness, false-detection rates, and confusion (see §3.2| §3.5| and 
§2.2.2| respectively). 

We populate the raw signal map for each cluster with simu- 
lated point sources of varying flux values at random positions and 
attempt to recover them with our standard source detection algo- 
rithm ( §2.2.2| ). The point sources are generated from the beam tem- 
plate and inserted one at a time to avoid significantly altering source 
density properties. The total number of sources injected per signal 
map is dependent on the real sources in that cluster; the upper flux 
limit is roughly a factor of two greater than the brightest cluster 
source. We inject 50 fake sources at each flux value, starting at 
0.5 mJy and stepping in 0.5 mJy bins. The positions are randomly 
chosen based on the noise mask criteria explained in §2.2.2| After 
implementing our cleaning and detection routine, we measure the 
output flux and position of the point source. If the source surpasses 
the 3 cr detection limit and is retrieved within 8 arcsec (roughly half 
the width of the SCUBA beam) of the input position, it is consid- 
ered recovered. 

We also repeat the simulations on maps where all the real 
sources have already been removed (i.e. any real source signal has 
been cleaned out). This is analogous to the criteria implemented by 
many groups (e.g. |Scott et al.|200'2l[2008j ) in which any fake source 



injected or retrieved within a certain distance of a real source is ex- 
cluded from further analysis. The flux boosting estimate (see §3.3| ) 
utilizes these simulations. 



3.2 Completeness 

The Monte Carlo simulations provide a means to estimate the com- 
pleteness individually for each cluster map, using the aforemen- 
tioned criteria for recovered sources. In Fig. [3] we plot the differen- 
tial completeness as a function of intrinsic source flux density for 
RCS 2319. Following ICbppin et al.| ([ 2006| , we fit to a curve of the 
form (5"*)/ {b + cS^), where S represents input flux. This empir- 
ically derived function provides a good fit to the calculated com- 
pleteness. For this specific cluster, we are ^50 per cent complete 
at 3 mJy and ^95 per cent complete by 9 mJy, which is consistent 
with the remaining six clusters. 



3.3 Flux Deboosting 

A well-established yet rather subtle bias, flux boosting describes 
the enhanced flux measurement for sources in low S/N (<10) maps 
that occurs when there is an increase in number counts of sources 
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Figure 3. Differential completeness versus intrinsic flux density for 3 a de- 
tections in the real signal map of RCS 2319. The error bars represent bino- 
mial uncertainties. The points are fit to a curve of the form (S'")/(6+cS'"), 
where S represents input flux. 



Figure 4. The flux boosting factor as a function of output flux (i.e. measured 
flux) for all cluster fields combined. All sources detected at a flux density 
where the boosting factor lies above the dotted line must be deboosted by 
that amount. 



with decreasing flux density ( |Hogg & Turner '' 1998). This pro- 
duces an asymmetric distribution of sources available for scattering 
into another flux domain. In other words, there is a greater proba- 
bility that intrinsically faint sources are detected at a higher flux 
than bright sources are deboosted to a lower flux, which is known 
as the Eddington effect. This measurement effect differs from a 
Malmquist bias, which arises in flux-limited samples due to greater 
detectibility of more luminous sources over a greater volume. A 
second systematic flux measurement error for low S/N sources re- 
sults from a tendency to be pulled to brighter fluxes due to random 
noise spikes in the map. This is a consequence of the method of 
recovery of sources since peaks define the position and flux of the 
source. 

To quantify this asymmetric flux error, many groups (e.g. 
ICoppin et aEHjOOSl [2006, Scott et al.|(2008] ) have implemented 
a Bayesian approach by calculating a posterior flux distribution 
(PFD). As an alternative to this analytical method, numerical simu- 
lations also provide a means to estimate the flux deboosting correc- 
tion (e.g. Sc ott et al.|20()2| ). A comparison of the output flux to the 
range of input fluxes that produced it should reveal a flux boosting 
factor that increases with decreasing measured flux density. 

Using the results from the injected point sources at each flux 
level within the uniform noise region, we bin all input fluxes (with 
1 mJy width) that yield a 1 mJy range of output fluxes for each clus- 
ter field separately. The centre of the output flux range is divided 
by the mode of input fluxes to calculate the flux boosting factor. 
Each cluster produces approximately the same results within un- 
certainties, therefore the same procedure is repeated combining the 
simulations from all clusters to increase number statistics. The final 
flux boosting factor as a function of output flux is shown in Fig. [4] 
This factor is then applied to the measured SMG fluxes within that 
output flux range, which are, consequently, deboosted. The boost- 
ing factor ranges from 1.7 in the lowest flux bins to 1.0 at higher 
fluxes. Sources detected above 6mJy are therefore left unaltered 
based on these simulations. The final deboosted 850 |^m flux den- 
sities are recorded in Table |2] 

3.4 Positional Uncertainty 

Lower S/N submillimetre sources are typically retrieved at larger 
positional offsets than more robust sources due to their suscepti- 



I2r^ — ' — ^ 
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Figure 5. The positional uncertainty of the SMG as a function of S/N based 
on our Monte Carlo simulations combined for all cluster fields. The points 
represent the 95th percentile of all positional offsets at that S/N. 

bility to being pulled off their intrinsic position by random noise 
spikes, or through confusion with brighter sources. Based on our 
simulations, we calculate the 95th percentiles for the positional off- 
sets after binning sources in S/N with 1 a bin widths. The resulting 
histograms of positional offsets for each cluster follow the same 
pattern, and we therefore combine all clusters to increase our num- 
ber statistics. We find maximum positional offsets ranging from 10 
to 2 arcsec for S/Ns between 3.5 and 7.5 a, respectively (see Fig.[5]), 
which are consistent with other submillimetre cluster surveys (e.g. 
|Webb et al.|2005l [Scott et al.|2008| ). These offsets are subsequently 
used as search radii for counterpart identification (described in 

3.5 False-Detection Rate 

The inclusion of the negative beams from the chop in the convolved 
image renders the instrumental noise symmetric about zero. How- 
ever, the confusion noise could potentially produce an asymmet- 
rical noise fluctuation distribution. We quantify this by populat- 
ing noise maps with sources distributed in flux according to the 
SHADES differential source counts ( [Coppin et al . "2006 ) fit to a 
Schechter function ( ||Schechter||1976| ), up to three times the RMS 
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Figure 6. The false-detection rate as a function of S/N in our combined 
cluster fields. The points are empirically fit to a power-law and rates are 
determined from this curve for spurious source corrections (e.g. in source 
counts analyses — see 



4 850 i^m SOURCE COUNTS 

A crucial goal of most submillimetre surveys involves an under- 
standing of the distribution of galaxies over a range of flux densi- 
ties in order to constrain evolutionary models. A plethora of studies 
have tackled this at 850 ^tm across a high dynamic range of fluxes 
(~ 0.1 mJy — 20mJy), utilizing massive cluster lenses to amplify 
the SMG population at the fa int end (e.g.|Smail et al.|1997 Cowie 



et al. '2002| |Smail et al.|2002l [Chapman et al. 2002 ; Knudsen et al. 
2006 ) and wide blank-field surveys to probe brighter SMGs (e.g. 
Scott et aLl [2002l [Borys et al.|[2503l [Webb et al.|[2003al |Coppin| 
et al.[2006l ). 

A comparison of the number of sources as a function of flux 
density (i.e. source counts) between cluster fields and unbiased re- 
gions of the sky (i.e. blank fields) can elucidate the influence of 
cluster environment on galaxy evolution through the verification 
of a statistical excess of SMGs in cluster fields. Our survey offers a 
further analysis within those cluster fields by examining differences 
in number counts between clusters with and without an abundance 
of optical arcs from strong gravitational lensing and probes higher- 
redshift clusters, exposing excess counts as due to either lensing or 
contamination by cluster galaxies. 



noise in the map. In other words, there should be no real, unblended 
sources with S/N > 3a. We then count both the positive peaks 
(>3cr) and negative peaks (<-3cr), average over 100 realizations, 
and find a symmetric distribution. For example, the RCS 2319 field 
yields 4.0 ±1.5 positive peaks, and 3.5 ±1.6 negative peaks. 

We therefore assume the noise is symmetric about zero, which 
allows for a direct method to quantify source reality. After remov- 
ing all the detected sources from the raw image, we invert the map. 
We can then estimate the false-detection rate by counting the num- 
ber of peaks ^ 3 a within the usable map area on the inverted im- 
age. To improve our statistics, we combine all the cluster fields in 
the calculation instead of considering them individually. 

Consistent with other groups (e.g. | Cowie et al.||2002| |Scott| 
|et al.|2002] ), we find a 3 a cut includes a high number of spurious 
sources, corresponding to a ^46 per cent false-detection rate. The 
situation remarkably improves for a threshold of 3.5 cr, with only 
three (of 19) spurious sources detected, resulting in a 16 per cent 
contamination. Fig. [6] plots the false-detection rate at average S/N 
values including sources down to 3.0 a. The points have binomial 
uncertainties and are empirically fit with a power-law of the form 
y — ax^ , where 941.5 and —6.7 are found for the amplitude and 
slope, respectively. The false-detection rate for a given S/N is de- 
termined from this curve in subsequent analysis. 

We can estimate our expected number of spurious sources by 
summing the false-detection rate equation over the total number 
of detected sources in our survey. We calculate that four of the 26 
sources in this survey should be false, given their range in S/N val- 
ues (see Table [2|. 

Despite the significant uncertainty in individual sources de- 
tected below 3.5 cr, we decide to include them in the source counts 
analysis (as do other groups — e.g. |Cowie et al.|2002{|Coppin et al. 
|2006| ) as a means for comparison to previous surveys. We correct 
the source counts for spurious sources (see and contend that 
the 3 cr threshold is reasonable for the counts analysis. We do note, 
however, that counterpart identification at this level can be difficult 
to accurately interpret, as individual sources carry low significance. 



4.1 Gravitational Lensing 

Strong gravitational lensing magnifies sources, and subsequently 
increases the observed flux. We therefore must assess the effect 
any flux magnification could have on the source counts. The cluster 
mass inferred from the richness estimate (see Table [T} can be used 
to estimate the magnification caused by the smooth cluster mass 
distribution. The amount of flux boosting depends on the source 
redshift, and as these are unknown, we consider values between 
z — 1 and z = 4. To compute the magnification, we assume that 
the density profile of the cluster is described by a Navarro, Frenk 
& White profile ( [Navarro et al.|1997| ), with a mass-concentration 
relation given by [Duffy et al. ( 2008| ). Given that the control and 
lensing samples span the same range of cluster mass, we treat both 
cases with only a smooth cluster component, and do not attempt to 
reproduce the strong lensing features described in §1.2| 

These calculations only consider the contribution of the 
smooth density profile, but locally cluster members themselves can 
contribute to the lensing signal as well. For SMGs in the outskirts 
of clusters, the source would have to be located very close to the 
lens to be magnified significantly. 

We find that the magnification is modest for most sources, but 
in a few cases the SMG is likely to be strongly lensed if it is indeed 
at high redshift and the adopted mass is representative. However, 
it is difficult to quantify without being able to model the cluster in 
more detail. 

Assuming z = 2.5, a typical SMG redshift, we calculate an 
average magnification of ^ 1.5, excluding the two SMGs that 
would be significantly lensed if confirmed to be at that redshift 
(SMMs 2319.3 and 1326.1). In particular, the magnification values 
for SMM 23 19.3 are unreliable as it is extremely close to the cluster 
centre (i.e. within a projected distance of 13 arcsec). This magnifi- 
cation value is fairly robust to 1 cr increases/decreases in the cluster 
masses, yielding average values of ^ 1.8 and 1.3, respectively, for 
di z = 2.5 SMG population. Even at the most extreme level, as- 
suming a source plane redshift of z = 4, we only find an average 
magnification of ^ 1.8. Due to the fairly low magnification values, 
we make no strong lensing correction in our source counts analysis 
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(discussed further in §4.2| ) and note that none of these systems have 
actually been confirmed as lensed objects. 

4.2 Differential Counts 

An accurate estimate of 850 |^m source counts relies on prop- 
erly correcting for several of the properties described above, such 
as incompleteness ( §3.2| ), flux-density boosting ( §3.3| l, and false- 
detection contamination ( §3.5| l. Following the procedure outlined 
in 'Scott et aLj ( |2002| ), we exploit the results of our Monte Carlo 
simulations to account for these effects. 

To begin, we devise flux density bins based on the values used 
in the completeness calculation of §3.2| (i.e. 0.5 mJy bins centred on 
integer and half-mJy fluxes). We then calculate raw source counts 
for each bin, separately for each cluster, by adding the number of 
objects that fall within that bin using their statistically deboosted 
flux density. To account for spuriously detected sources, each raw 
count is weighted by the probability that it is false, which is based 
on the S/N of the source. The probability is estimated from the 
fit to the curve in Fig. |6] This weighting scheme de-emphasizes 
counts from lower S/N sources which are more likely spurious. 
Due to the observed incompleteness of each map, which increases 
with decreasing flux, we adjust the weighted counts by dividing 
them by the best-fitting completeness value at that flux determined 
from a curve analogous to that in Fig. [3] (each cluster field has a 
slightly different completeness fit). Bins with zero initial counts 
are left unaltered. These corrected values are combined, preserv- 
ing bins, for each subsample separately (i.e. lensing and control 
groups). Based on the distribution of flux bins, we create larger 
bins to ensure all differential bins have a value greater than zero for 
display purposes. This results in four new bins for each subsample, 
which are presented in Table [3] The highest flux bin in the control 
group is a result of one extremely bright source in RCS 2318, de- 
tected at ^^15 mJy. After adding the values from all the old bins that 
fit within the new bin boundaries, we divide each bin by its width 
for normalization purposes. The final differential counts, shown in 
Fig. |7] are simply these values divided by the total area in square 
degrees covered by that subsample (see §2.2.2| for a description of 
how usable map areas are determined and Table [T] for a list of ar- 
eas). The relative uncertainties are calculated from standard Pois- 
son errors on the raw counts scaled by the bin size and map area. 
The final counts and errors are tabulated in Table [3] 

In Fig. |7] we also plot the differential counts obtained from 
the SHADES blank-field survey, fit with the empirically motivated 
Schechter function ( |Schechter|1976| ) 



Table 3. Differential 850 |xm source counts for the control and lensing sam- 
ples. 



dN 

Is 



(1) 



The SHADES team has found the best-fitting form to the function 
using values of 3.3, 1599, and -2.0 for S\ N\ and a, respectively. 
The differential counts of our control sample agree quite well with 
the SHADES results; the lensing sample tends toward a slight ex- 
cess in source counts compared with both the control sample and 
SHADES survey counts, hinting at an inherent difference between 
the cluster environments of the two subsamples in our survey. We 
note that demagnifying the fluxes of the lensing cluster sample by 
a minimum of 1.5 does bring the differential counts in line with 
SHADES. However, if we assume the field is lensed, we also need 
to decrease the total survey area, since massive foreground lenses 
cause the background area to look expanded ( |Broadhurstl| 1 995| ) . 
This brings the counts back up, and thus we emphasize that this is 
only a lower limit on the required magnification. Since the control 
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Figure 7. Differential 850 \im source counts. Our results are denoted by the 
larger symbols, where the circles represent counts from the lensing clusters 
and the diamonds represent counts for the control group. We also include 
differential counts from the largest blank-field survey — SHADES (crosses) 
jCoppin et al.|2006} — and a best-fitting Schechter function. Our larger error 
bars are attributable to small number statistics. 



sample covers the same range of magnifications, the same adjust- 
ments would be made for both samples, and the counts discrepancy 
between the two would remain. As none of these galaxies have been 
confirmed as lensed systems, it is difficult to fully assess the mag- 
nification correction. 



4.3 Cumulative Counts 

The same corrections (i.e. flux deboosting, false-detection rates, 
and completeness adjustments) are applied to properly compute 
source counts at each 0.5 mJy flux density bin. The cumulative 
counts are recorded at each discrete flux density as the total num- 
ber of sources per square degree at that flux and higher, meaning 
each bin is correlated with all higher flux density bins. The control 
sample ranges from 1.0 to 5.0 mJy, while the lensing group encom- 
passes fluxes between 2.0 and 8.0 mJy. The resulting counts are 
given in Table|4]and shown in Fig.|8]along with various blank-field 
and cluster surveys. 

We test the validity of our source count results by randomly 
eliminating one cluster from each sample and re-calculating the cu- 
mulative counts. Separately, we also transfer one control cluster to 



© 2011 RAS, MNRAS 000,[l]j28] 



SMG source counts in high-z galaxy clusters 1 1 



Table 4. Integral 850 \im sources counts for the control and lensing samples. 



Integral Counts 
N{> S) (deg-2) 
Flux Density (mJy) Lensing Sample Control Sample 
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7.0 
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the lensing sample to ensure the discrepancy in source counts be- 
tween the two samples does not flip. The cumulative counts prove 
to be robust to all test cases as the lensing sample counts remain 
above the control sample. We also attempt to account for any hu- 
man biases (mentioned in that may enhance the source counts 
in the lensing sample. We note that the sample was selected in 
an unbiased way, and no clusters have been left out of the analy- 
sis. As evident in Figures [T] and |2] the lensing clusters have larger 
map areas as the result of multiple pointings. These three clus- 
ters were all mapped in classical observing time, and therefore 
observations were designed in real-time to maximize the S/N of 
detected sources; the control clusters, on the other hand, were all 
observed later in service time. To determine the potential effect this 
bias might have on our results, we compute the cumulative source 
counts within the same projected radius from the original point- 
ing for each cluster. We adopt a radius of 80 arcsec to ensure the 
search area is within the uniform noise for each field. The resulting 
counts are plotted alongside the original counts in Fig. [8] We find 
our results are robust to any potential bias resulting from multiple 
pointings, and we therefore continue our analysis with the original 
source counts from the larger map areas. 

Although just below the SHADES results, the control sam- 
ple is consistent with the counts from the Large Apex Bolometer 
Camera (LABOCA) Extended Chandra Deep Field South (ECDFS) 
Submillimetre Survey (LESS); the ECDFS has a deficit of bright 
SMGs and is known to be an underabundant ULIRG field ( WeiB 
|et al.||2009| ). We note that we have plotted the LESS counts esti- 
mated directly from the catalogue, rather than those estimated from 
the P{D) analysis, although both methods show a deficit of counts 
compared to other deep fields. The lower counts for both LESS 
(with the direct source count method) and our control group could 
be the result of faint multiple sources in the beam identified as a 
single brighter source, which would lead to an underestimate of 
fainter sources ( WeiB et al. 2009). 

With the exception of the highest flux bin, the lensing sample 
shows an excess of counts compared to blank-field surveys, con- 
sistent with other lensing cluster surveys (e.g. |Cowie et al.|2002[ 
[Smail et "aL) l2002). This is a ^ 2.7 a overabundance relative to 
the control sample (at the lowest flux point, which accounts for all 
sources with S ^ 2 mJy), suggesting that enhanced source counts 
are unique to strong lensing cluster fields. We limit our quantitative 
comparison to the two samples presented in this work, as we can 
be assured that the counts were calculated identically, and do not 
suffer from any systematic differences. 



5 COUNTERPART IDENTIFICATION 

A wealth of ancillary data to the submillimetre observations exists 
at mid-infrared (mid-IR), radio, optical, and x-ray wavelengths (see 
§2.3| ), allowing for the identification of counterparts to the SMGs 
which improves positional accuracy and helps constrain redshifts. 
A possible complication with regards to identifying counterpart 
emission with ancillary observations involves the detection limit at 
various wavelengths. To ensure we are sensitive to detecting coun- 
terparts, we plot the 3 a depths for the MIPS, IRAC, and radio fields 
in Fig. |9] overlaid with two spectral energy distribution templates 
for Arp 220, a typical star-forming ULIRG at low redshift, and for 
M82, a star-forming galaxy with irregular morphology. Each SED 
has been converted to the highest redshift for the entire sample, 
z ^ 1, and normalized to the SCUBA detection limit at 850 \im 
in our cluster fields, 3cr 3mJy. The detection limits for 3.6, 
4.5, 5.8, 8.0, and 24 |^m along with 1.4 GHz are below the SEDs, 
thereby allowing sensitive counterpart detection for cluster mem- 
bers. A non-detection may by indicative of either a high-redshift 
background source or a spurious SMG, particularly for low S/N 
cases. We reiterate that although we expect four spurious sources 
given our false-detection rate equation (see §3.5| l, we search for 
counterparts in the entire source catalogue since we do not know 
which of these are in fact spurious. However, we limit our potential 
cluster members (see §6.1| l to sources with S/N>3.5. 

Given the fairly homogeneous data coverage with MIPS 
(imaging over six of the seven cluster fields), its higher positional 
accuracy (5.7 arcsec point spread function) compared to SCUBA, 
and the low density of mid-IR sources compared to optical, we take 
the MIPS imaging to be our primary dataset for counterpart identi- 
fication. The radio data, which covers four cluster fields, is used as 
a secondary indicator. We classify our identifications as secure (ro- 
bust 24 i^m detection with supporting radio and/or IRAC emission), 
ambiguous (the 24 \im emission present suffers from high confu- 
sion or there are two equally robust identifications), or tentative 
(there are a few possible candidates, but none are overwhelmingly 
convincing). The optical and X-ray data are only used to help es- 
timate redshifts and confirm the astrometry; they do not contribute 
to the strength of the classification. 



5.1 MIPS Counterparts 

The utility of mid-IR follow-up observations in the identification of 
SMG counterparts has been exploited by many groups (e.g. Ivison] 
let aL|2004[[Pope et al.|2006l [Biggs et aL|2010| ). Like the submil- 
limetre regime, mid-IR imaging probes dusty star-formation and 
AGN activity by sampling thermally emitting dust grains. Obser- 
vations of high-redshift clusters at 24 \im also provide a rough es- 
timate of the star formation rate through an extension of an empiri- 
cal correlation between it and the infrared luminosity of a starburst 
galaxy CKe nnicutt|1998bl|Bell|2003l|Calzetti et al.|2007[|Kennicutt| 
et al. 2009). 

To locate counterpart 24 ^im emission, we employ a signal- 
to-noise dependent search radius using our positional uncertainty 
simulations described in §3.4| which represent the 95th percentile 
of SMG offsets. We limit the maximum radius to 8.0 arcsec (since 
this defines the source recovery cut-off in our simulations from 
§3.1| l and the minimum radius to 2 arcsec, where the positional off- 
set tapers off in Fig. [5] We also introduce an additional 2.0 arcsec 
to account for uncertainties in the MIPS source position, resulting 
in a final search range up to 4.0 arcsec for the highest S/N SMGs 
and 10.0 arcsec around the lowest S/N SMGs. 
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Figure 8. Cumulative 850 fxm source counts for the lensing clusters (large red circles) and the control clusters (large blue diamonds). Also shown are counts 
determined from other cluster surveys jSmail et al. 2002 Cowie et al. 2002) and blank-field surveys ( Scott et al. 2002 Webb et al. 2003a Coppin et al. 2006] 
|WeiB et al.|2009j . We have scaled the 870 |j.m LESS counts to 850 |j.m using a submillimetre spectral index of 2.7. Note the tentative excess of source counts 
for the two cluster surveys mentioned, along with the counts derived from the strong lensing sample. 
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Figure 9. SED templates for Arp 220 (solid line) and M82 (dotted line) red- 
shifted to 2; = 1.0, the highest redshift of our submillimetre cluster sample. 
The templates have been normalized in flux density to 3.1 mJy (roughly 
the 3 cr detection limit in our SCUBA maps) at 850 fxm. The 3 a detection 
limits for the MIPS, IRAC, and radio data are also shown. Note that these 
local galaxy SED templates are not representative for high-redshift SMGs 
at rest-frame optical and UV wavelengths (i.e. below ~1 fxm). 



We find single MIPS counterparts for 13 (^ 57 per cent) of 
the 23 SMGs with 24 \im coverage. Three SMGs lack any 24 \im 
emission, and seven have multiple (two or three) MIPS identifica- 
tions within the search radius (see §5.1.1| l. The complete counter- 
part identification results are shown in Table |5] 



5.1.1 Multiple MIPS Counterparts 

In the case where multiple 24 \xm sources are found within the 
search radius (35 per cent of the 20 sources with mid-IR emission), 
we investigate whether one source in particular dominates the in- 
frared flux, and if so, take this to be the more likely counterpart. 
[Pope et aL] ( [200 6') conclude that typically only one source signif- 
icantly contributes to the submillimetre flux, while the secondary 
sources within the search area are incapable of producing the flux 
due to faintness. 

For the purposes of this work, we take the brightest MIPS 
source within the SCUBA beam to be the dominant emitter, but 
note that this could potentially bias us towards foreground galax- 
ies. In all cases, the brightest MIPS source is coincident with emis- 
sion at another wavelength if there exists any within the SCUBA 
beam, and therefore provides evidence in favour of our counterpart 
assumption. There is one source in RCS 1419 that has two MIPS 
sources of roughly equal fluxes (and no emission at any other wave- 
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Figure 10. A plot of the likelihood of finding one or more source within an 
aperture of the given radius. We linearly interpolate between these values to 
assign likelihoods to MIPS counterparts based on its offset from the SMG. 

length) for which we take the closer of the two (for both 850 and 
450 \im) as the correct identification. 

5.1.2 Assigning MIPS Likelihoods 

To quantify the probability that the MIPS detection is associated 
with the 850 \im emission, we perform Monte Carlo simulations, 
placing apertures of various radii on each 24 \im cluster field and 
counting the number of sources detected within that area. The im- 
portance of using our own data for the simulations, rather than a 
pure statistical method such as the P-statistic (e.g. Downes et al.| 
|1986| ) which requires an estimate of the source density, derives 
from the prior expectation of an overdensity in our cluster fields. 
We drop 1000 apertures at each radius from 1.0 to 10.0 arcsec in 
steps of 1.0 arcsec on the MIPS images (see Fig. [T0|. A likeli- 
hood for each radius is assigned based on the number of times we 
detect one or more MIPS sources within that aperture. Using the 
positional offset between the SCUBA source and the most robust 
MIPS counterpart (in the case of multiple 24 |^m sources within the 
search radius), we interpolate between the adjacent aperture radii 
(see Fig. [To]) to determine an individual likelihood that each MIPS 
counterpart is not just a chance alignment. 

To investigate the importance of random associations that con- 
taminate our identifications, we compare the positional offsets be- 
tween secure counterparts and their associated SMG to our Monte 
Carlo simulations of submillimetre positional uncertainties ( §3.4| ) 
combined for all clusters (see Fig.|ll|). The simulations have been 
weighted proportionally by S/N based on the percentage of SMGs 
with MIPS counterparts in each S/N bin (with 1 a width). This 
places more emphasis on lower S/N (i.e. 3-4 a) positional un- 
certainty simulations as 60 per cent of our submillimetre detec- 
tions fall within this range. We also include the theoretical distri- 
bution of radial offsets for flux-boosted SMGs from ilvison et al.l 
( [2007 Assuming a Gaussian distribution of errors in right as- 
cension (RA) and declination, this is given by re~^ /^'^ , where 
a = 0.66'[S/N^pp - (2/3 + 4)]-^/^ and the S/Napp has not been 
corrected for flux boosting. 

The simulated and theoretical distributions agree quite well 
and both peak around 2.5 arcsec, while our broader distribution of 
MIPS counterpart offsets peaks further out at ^ 6 arcsec. This sug- 
gests that some of our counterparts are falsely identified, primarily 
those with offsets beyond 6 arcsec, where 50 per cent of the his- 



Figure 11. The distribution of positional offsets for the brightest MIPS 
counterparts in our sample (dashed histogram) compared to Monte Carlo 
simulations of the SMG positional uncertainties (solid histogram). We also 
overplot the theoretical distribution of radial offsets for flux-boosted SMGs 
(dotted curve), as defined in Ivison et al. 2007 using the mean S/N of 
the SMG sample. The dash-dotted histogram at 11-12 arcsec denotes the 
fraction of unrecovered submillimetre sources (recovered at offsets > 8 
arcsec) in our simulations, meaning 15 per cent of our sources should not 
have counterpart identifications. Note that this distribution is for all 850 |j.m 
sources with S/N ^ 3.0. 



togram of MIPS counterparts lies above the simulated histogram 
(including unrecovered sources). Of the ten MIPS sources in this 
high offset region, five are securely identified (see ^for a defini- 
tion) as the correct SMG counterpart, while the remaining 50 per 
cent are only tentative or ambiguous counterparts (see Appendix 
A for details on classifications of counterparts). Our classifications 
are therefore consistent with the distributions at large positional off- 
sets, with tentative counterparts more likely being wrongly iden- 
tified. Alternatively, these sources with large MIPS offsets could 
also represent blends of multiple submillimetre sources, or spurious 
sources (see §3.5| ). In fact, four of the sources with MIPS offsets 
above 6 arcsec have S/N<3.5, increasing the likelihood that they 
are either spurious or blended sources. Moreover, the fraction of 
unrecovered submillimetre sources in our simulations (dot-dashed 
histogram) provides a handle on how many real sources should re- 
main unidentified, roughly 15 per cent, which is congruous with 
our findings of three SMGs that lack MIPS counterpart emission 
within the search radius. 



5.2 Radio Counterparts 

The powerful benefits associated with complementary radio data 
at 1.4 GHz (21 cm) with the VLA are primarily two-fold. First, the 
sub-arcsecond spatial resolution attainable with radio interferome- 
try provides improved astrometric precision compared to SCUBA 
and MIPS positions. Second, radio sources are fairly sparse in the 
sky ( |Richards|2000| ), making coincident detections of 850 ^im and 
1.4 GHz emission highly unlikely. Furthermore, radio emission is 
sensitive to synchrotron radiation from relativistic electrons asso- 
ciated with supernovae, and thus traces recent star formation, re- 
sulting in a correlation between radio and submillimetre emission. 
Any radio emission detected within the 850 |^m search area, there- 
fore, offers likely counterpart identifications to the SCUBA sources 
( Ivison et al. 2000). In principle, the simulations for the MIPS like- 
lihoods ( §5.1.2| ) could be applied to the radio maps as well, but be- 
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cause the number of radio sources is even sparser, we automatically 
assign the 1 .4 GHz emission as the correct identification. Note that 
one exception to this occurs in source SMM 2319.3, in which two 
radio galaxies are found; further discussion is provided in §6.3| 

Radio counterparts are searched for within a maximum of 
Sarcsec radius around the SMG, depending on the S/N of the 
source. No additional margin is given due to the high astromet- 
ric precision of the radio data. For the four clusters (17 SMGs in 
total) with complementary 1 .4 GHz imaging, we detect six radio 
sources within the search radii of five SMGs, although one source 
is suspected to be unassociated with the SMG (discussed in detail 
in §6.1.1| l In all radio-detected cases, a secure 24 \im counterpart 
has already been identified, and the MIPS and radio positions are 
within 3 arcsec of each other. 



|2000 | have exploited hard (i.e. 2-8 keV) X-ray observations to lo- 
cate any AGN associated with submillimetre emission which might 
contribute to the extreme luminosities produced by these sources. 

X-ray counterparts must adhere to the same distance criteria 
as radio counterparts — within a maximum of 8 arcsec of the SMG 
(depending on the S/N) and 3 arcsec of a secure MIPS counterpart. 
Only five X-ray sources are found within the search radii, leaving 
16 SMGs without any X-ray association. 



5.3 IRAC Counterparts 



Many groups (e.g. Smail et al. 1999 , Ivison et al. 2002 , Webb et al. 



2003b[[Ashby et al. 2006, Pope et al. 2006, Yun et al. 2008, Hain 
line et al.|2009 JWardlow et al.|2010| ) have investigated the efficacy 



of infrared colours of SMGs as a rough proxy for redshift as well 
as a counterpart selection technique. As these colour criteria can be 
biased towards higher redshift SMGs and in some cases, are per- 
haps too general (as found by ,Hainline et al.|[2009] ), we employ 
them merely as a means to secure one identification over another in 
the case of multiple MIPS detections within the search radius. 

Previous studies (e.g. |Yun et "aT] |200 8) have used a 
S's.o/xmZ&.s/xm vcrsus S5.8 ^^m / S3.6 colour-colour diagram 
(see Fig.[T2| to identify mid-IR counterparts to SMGs based on the 
claim that their colours occupy a different region of the plot than 
those from stellar photospheric emission (the cloud in the blue re- 
gion) and of low- 2; late-type galaxies (the vertical branch emerging 
from the cloud). These colour criteria derive from the region popu- 
lated by AGN identified byl Lacy et ahlpOOS] !. In Fig. [12] we show 
the infrared colours for both secure and tentative IRAC counter- 
parts to the SMGs against a backdrop of galaxies in the six cluster 
fields. Overlaid are SED tracks for representative galaxies, includ- 
ing the ULIRG Arp 220, an irregular starburst galaxy, an elliptical 
galaxy, NGC 1068 (an AGN), and a pure power-law galaxy. Of our 
secure counterparts, 11 of 12 (92 per cent) fall, within the errors, 
inside the marked SMG domain. The remaining secure counterpart 
falls along the low-^ late- type sequence and has only a single MIPS 
detection in the entire search radius; it is therefore still identified as 
a secure detection. Thus, the IRAC observations support all of our 
previous conclusions based solely on the MIPS and radio data. 

We also plot a colour-magnitude diagram (Fig. [T2| for the 
same counterparts described above using the criteria suggested by 
[Pope et ar] ( |2006| ) and |Hainline et al.| ( [2009| ). This domain is used to 
select z > 1.5 SMGs with a 90 per cent success rate and 6 per cent 
contamination rate by z < 1.5 galaxies (based on the more relaxed 
[Hainline et al.|2009l criteria). This diagram also serves as a selec- 
tion cut for SMG counterparts as there is a low probability of ran- 
domly displaying these colours. Regarding our secure counterparts, 
75 per cent are within error (la error bars on infrared colours) of 
the z > 1.5 region, suggesting the majority of the SMGs in the 
sample lie at higher redshifts and are not cluster members. 



5.4 X-ray Counterparts 

The X-ray regime provides a window into obscured AGN and star 
formation; many groups (e.g. [Almaini et al.||199^ [Fabian et al.| 
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Figure 12. [LeftJS'g.o /^m/*S'4.5 versus 5*5.8 ^m/Ss.Q ^rn colour-colour diagram for both secure (large filled circles) and tentative/ambiguous (large empty 
circles) IRAC counterparts to the SMGs. The filled stars denote likely cluster members (see §6.1.1} . Small black dots correspond to field galaxies in the six 
clusters detected above 5 cr in all four IRAC channels. SED tracks for representative galaxies are overlaid from z = Otoz = 3.1, with the triangle indicating 
low redshift. Also shown is a pure power-law galaxy, S (x. , with a spectral slope ranging from —0.5 to —3.0. The cluster redshifts of RCS 0224 (0.113) 
and 2319 (0.9024) are displayed along the tracks as red and cyan squares, respectively. The dotted gray line represents the boundary for separation of z > 1 
SMG counterparts from field galaxies as suggested by |Yun et'dr| 2008 [Right] A colour-magnitude diagram for the secure and tentative IRAC counterparts. 
The symbols are equivalent to those in the colour-colour diagram to the left. The solid black line denotes the region proposed by |Pope et al.|2006) for SMG 
counterpart selection, while the dotted gray line traces the more relaxed criteria given by |Hainline et al.|2009| 
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6 DISCUSSION 

6.1 Tentative Cluster Members 

After careful inspection of all counterpart identifications, we infer 
that four (SMMs 0224.2, 0224.3, 2152.2, and 2319.3) of the 26 
SMGs are potential cluster members based on their photometric or 
spectroscopic redshift. We note that this conclusion is biased to- 
wards cluster fields with radio and/or multi-band optical data, as 
these wavelength regimes provide positional accuracy and redshift 
estimates. With more consistent wavelength coverage for all clus- 
ters, the fraction of possible members could be higher. 

For the remainder of this paper, however, we choose to limit 
our discussion to the three most likely cluster members and elimi- 
nate one candidate (SMM 2152.2) based on its low (< 3.5) S/N at 
850 i^m. As determined in §3.5| a detection threshold of 3.0 a car- 
ries a high false-detection rate and therefore these sources are not 
reliable as individual objects. Three remaining candidates (SMMs 
0224.2, 0224.3, and 2319.3) form the basis of our cluster member 
sample and are discussed in detail here. Multi- wavelength postage 
stamps of the counterpart emission for these three SMGs are shown 
in Fig. [13] A discussion of counterparts for the remaining SMGs is 
deferred to the appendices. 



6.1.1 Notes on Individual Cluster Members 

SMM 0224.2. This SMG has a secure identification based on 
a robust 24 ^tm detection. Coincident radio. X-ray, IRAC, and 
optical emission all fall within 0.5 arcsec of the MIPS position and 
greatly strengthen the security of this identification (see Fig. [13]). 
The absence of 450 ^im emission is likely due to the high noise 
at the edges of the 450 ^tm FOV where this SMG is located. The 
slight offset between the optical and radio centroids is probably 
due to the low S/N (3.8 a) of the radio detection. The infrared 
emission is significantly detected in all four IRAC channels and 
falls within the SMG counterpart region of the S5.8 fim/ S3.6 fim 
versus S^.g^m plot, making it a likely SMG counterpart (see 
[Pope et al.|[200gl [Hainline et~ari 12009] ). Although this region 
of the colour-magnitude plot typically accounts for z > 1.5 
SMGs, the S5.8 fim/ S3.6 nm colour could be reddened due to AGN 
contamination. Indeed, the source is along the power-law AGN 
region in Sg.o fim/ S4.5 versus S5.8 fim/ S3.6 fim colour-colour 
plot, and has MIPS -IRAC colours indicative of a starburst (SB) + 
AGN component (see §6.1. 2| ). The photometric redshift places the 
object at 2; = 0.677 ± 0.074, consistent with the cluster redshift at 
the 1.3 cr level (see also |Webb et aL|2005| ). 

SMM 0224.3. Although two sources of 24 \im emission are 
found within the 10 arcsec search radius, their flux densities differ 
by a factor of 2.5. At 404 |aJy, the brighter (but further) source co- 
incides with X-ray emission and an IRAC source robustly detected 
in all four channels, further securing the correct identification. 
The infrared colours of Sg.o ij.m/ S4.5 ij.m and S's.s/xm/^'s.e/xm fall 
within the likely SMG counterpart region of the colour-colour plot 
as denoted by |Yun et al.H2008] ) and are marginally consistent with 
the cluster redshift on the irregular starbursting galaxy track. It also 
falls within the appropriate SMG domain proposed by |Hainline| 
|et al.| ( [2009] ) on the colour-magnitude diagram. We note there is no 
1.4 GHz detection for this SMG; weak radio emission is found to 
the southeast, but with too large an offset from the MIPS position 
to be associated. The lack of 450 \im emission is a result of the 
noisy edges on the map. The photometric redshift for a coincident 



optical source places this galaxy at the redshift of the cluster (at 
the 1.2 cr level), with z = 0.688 ± 0.071. 

SMM 2319.3. SMM 2319.3, the most intriguing SMG of the 
sample, represents the only field where two radio sources have 
been detected. This is not necessarily indicative of false radio 
detections as 10 per cent of SMGs are expected to contain two 
radio sources within the search radius, both of which are related 
to the 850 |j,m em ission pvison et al.l2002[ [Chapman et al. 12005] 
[Pope et al.|[2006] [Ivison et al.pOOvf . Furthermore, this SMG is 
positioned very close to the cluster centre (^-^ 10 arcsec) where the 
density of galaxies is high. 

Both radio sources are extremely bright at 1 159 and 5656 pJy. 
Closer investigation reveals a possible head-tail morphology for the 
brighter object; it is an extended radio source with a single bright 
optical source at one end. It is most likely the result of radio jets 
bending due to the bulk flows in the intracluster medium, and pos- 
sibly suggestive of merger-induced activity. It has a spectroscopic 
redshift consistent with the cluster redshift {z — 0.90490). How- 
ever, the spectrum shape is indicative of an old elliptical galaxy, 
which is not typical for SMGs. The fainter radio source is also 
a spectroscopically confirmed member {z = 0.90126) and most 
likely associated with the brightest cluster galaxy (BCG). Spectro- 
scopic observations were obtained with the Visible Multi-Object 
Sectrograph (VIMOS; [Le"Fevre et al.|2003 ) on the 8.2-m VLT UT3 
telescope, and the Inamori Magellan Areal Camera and Spectro- 
graph (IMACS; |Dressler et al.|2006| ) on the 6.5-m Baade telescope, 
for the head- tail source and BCG, respectively. 

The radio sources lie on opposite sides of a 24 \im source, 
within 1.9 and 2.7 arcsec, for the fainter and brighter source 
respectively. Given the above description on the spectrum shapes 
and positional offsets, the MIPS emission is likely associated with 
the fainter radio source. The IRAC emission is highly confused, 
complicating the identification process. X-ray emission is detected 
as well, but at a distance of 3.0 arcsec from the fainter radio source; 
it is more closely aligned with the MIPS emission. The slightly 
lower S/N of the X-ray source (3.9 a) might explain the larger 
offset between it and the more accurate radio position. Another 
possible interpretation could involve an ongoing merger event 
between the MlPS/X-ray and radio sources, as there are only 
^15kpc between the MIPS and radio centroids. Based on the 
proximity to the cluster centre, the submillimetre emission could 
also be a background lensed source; however we have no evidence 
to substantiate this possibility. The SMG search radius also has 
three optical sources with photometric redshifts consistent with 
the cluster. Although the crowded field complicates the process 
of identifying coincident emission, it seems clear that this is a 
submillimetre detection of the cluster centre, and therefore a likely 
member. 

Thus far, only a handful of SMG studies have detected submil- 
limetre cluster galaxies, some of which where discovered serendip- 
itously as a result of utilizing massive clusters as gravitational tele- 
scopes to probe the background population (e.g. jSmail et al.|2002] 
[Cowie et al.|2002| ). Other studies have indeed focused on the cluster 
galaxies, but either rely on unconfirmed members (e.g. Best"20021) 
or study lower redshift clusters than presented here, for example 
two clusters at z = 0.23 and z — 0.25 that each contain a cen- 
tral galaxy with submillimetre emission ( [Edge et al.|1999] ), and a 
rich cluster at z = 0.54 with two likely 1.1mm detected cluster 
members ( Wardlow et al.|20l"m . Our detection of three likely sub- 
millimetre cluster members at the higher redshifts of ^ = 0.773 
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Figure 13. 24 x 24 arcsec^ postage stamps for the three cluster members: SMGs 0224.2 (top), 0224.3 (middle), and 2319.3 (bottom). The grey-scale images 
are optical z' band, IRAC 3.6 |Lim, and MIPS 24 |Lim, from left to right, where north points up and east is to the right. Each image is centred on the SMG 
position, denoted with a large blue cross. The size of the cross represents the search radius for MIPS counterparts. The rest of the symbols are as follows: 
a yellow square is a MIPS source, magenta triangles signify IRAC detections, the red contours show radio emission, an orange cross is 450 |Lim emission 
(although none is present in the candidate members), purple diamonds show optical sources with photometric redshifts, and a green 'x' denotes an X-ray 
detection. Note that the radio contours cover the entire field, but the symbols are only shown for sources within the search radius. The radio contours start at 
45 |Ljy and increase in 10 fjJy steps up to 105 fjJy, after which they increase by 200 [jJy. 



and z = 0.9024 not only has implications for galaxy cluster evo- 
lution, but also serves as an important probe of some of the most 
intense starbursts in clusters. 



6.1.2 AGN Contamination 

Infrared emission is a well-established tracer of star formation ac- 
tivity, in particular for very luminous and dusty galaxies. An excep- 
tion to this involves the contamination to the infrared luminosity 
due to AGN activity. 

In order to isolate sources contaminated by reprocessed AGN 
emission, we employ various infrared colour selection techniques 
(e.g. [Ivison et al.l[2004l |Lacy et al.||2004| |Ashby et al.1[2506l 
|Yun et al.||2008[ iHainline et al.| |2009 ). For instance, the left 
panel in Fig. [12] displays the redshift track of an AGN in colour- 



colour space along with that of a pure power-law galaxy, S oc 
jy"^. Of the three cluster members, SMM 0224.2 is the most 
likely AGN candidate, inhabiting the power-law AGN region, with 
log(5'8.o/xm/S'4.5/xm) = 0.21, wMch is consistcnt with AGN 
colours and redder than the majority of SMGs which typically have 



a ratio < 0.2 (see |Yun et al. 2008 ). Furthermore, according to|Ivi-j 
[son et al.| ( |2004| ) and Hainlin e et al.| ( |2009| ), the colour ratios of 
S24f,m/S8f,m (3.3) vcrsus Sg .0 f^m / S4.5 (1-6) for RCS 0224.2 
are indicative of a starburst (SB) + AGN component. 

In contrast, SMM 0224.3 shows no evidence for AGN con- 
tamination. A comparison of the MIPS and IRAC fluxes yields a 
starbursting system with no AGN component (i.e. S24 ij,m/ Sg nm = 
16.2 versus Sg.o f^m/ — 0.8; see Ivison et al.|2004l l. In the 
IRAC colour-colour plot, this SMG lies directly between the two 
starbursting tracks. SMM 2319.3 is too confused in IRAC to accu- 
rately interpret its infrared colours. 
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Our AGN analysis is further supported by empirically fitting 
the IRAC SEDs of the two cluster members with good infrared 
detections to a power law (e.g. |Ashby et"ar]|2006[ [Barmby et al.| 
12006] ). The sole AGN candidate, SMM 0224.2, can be fit to a red 
power law with a ~ —1.2, and monotonically increases in flux 
with wavelength through the IRAC bands; these features are both 
indicative of AGN-dominated emission. SMM 0224.3 has a blue 
power law (a ^ +0.4), suggesting stellar-dominated emission, and 
SMM 2319.3 is too blended in the IRAC channels to analyse. 



6.1.3 Cluster Member Properties 

We estimate the infrared luminosities (over 8-1000 ^im) of our 
cluster members following the prescription outlined in |Elbaz et al.| 
( [20Q2l l and [Sanders & Mirabel] ( p^96l ). Using a galaxy SED tem- 
plate for Arp 220, we shift the distribution to the respective red- 
shifts of the host clusters, normalize to the 850 \im flux of each 
SMG, and convert to a total infrared luminosity. This yields Lir = 
6.3 X 10^^ Lq, 5.6 X 10^^ Lq, and 3.9 x 10^^ Lq for SMMs 

0224.2, 0224.3, and 2319.3, respectively, which are consistent with 
ULIRG level luminosities, albeit on the higher end of the distribu- 
tion. 

Luminous galaxies (such as SMGs) powered by dense star- 
bursts output the bulk of their luminosity in the infrared as young 
stars heat the dust, which thermally re-emits the radiation. The in- 
frared luminosity, therefore, traces young stellar populations and 
provides an estimate of the star formation rate (SFR). Using the 
relation derived in |Kennicutt| ( |I998a| ) and assuming pure star- 
bursting systems, we calculate star formation rates corresponding 
to 1082, 968, and 665M0yr"^ for SMM 0224.2, 0224.3, and 

2319.3, respectively. This is in good agreement with rates derived 
in larger surveys with spectroscopically (or photometrically) iden- 
tified SMGs (e.g. [Chapman et"ari[2005| [Pope et al.[[2006| , where 
they calculate a median SFR of IIOOMq yr~^. We note that our 
estimates are highly sensitive to the dust temperature (Webb et al.' 
[2005 ), and could also suffer from AGN contamination, in particular 
SMM 0224.2, which would bring down our estimates. 

The intense star formation in these likely cluster members 
prompts the broader question: what is the trigger for the creation of 
these luminous systems? We can attempt to draw clues from their 
radial distribution. For instance, the potential member within RCS 
2319 is located in the cluster core; the MIPS emission is ^ 8kpc 
from the brightest cluster galaxy (BCG; see Fig.[TJ. Such high ac- 
tivity in the proximity of the core is indicative of a major merger- 
induced episode of star formation as a bound pair falls through the 
cluster centre and is tidally torqued by the potential (e.g. [Oemler[ 
[etal.|2009] ). Moreover, recent evidence suggests that the build-up of 
BCG mass occurs via major mergers f Tran et al.|2008 i Either sce- 
nario could explain the extreme star formation in the core, though 
the high-density of sources complicates identification of counter- 
part emission. 

In contrast, the two members of RCS 0224 are both further 
out at ^ 95 arcsec from the BCG, corresponding to ^ 700 kpc 
at z = 0.773, and therefore could belong to the infalling popula- 
tion. Indeed, galaxy infall events have been observed to increase 
at higher redshifts (e.g. [Filing son et al.|200l1[Loh et al.|2008l ). In 
this scenario, the increased star formation of these accreted systems 
might merely reflect the increased activity in the field, as suggested 
by numerous studies (e.g. Tran et al.|2005[[Saintonge et al.|2008| ). 



6.2 Redshift Estimates 

For the remaining 23 SMGs in the sample, the counterpart identifi- 
cation process revealed only two secure identifications with either 
an optical spectroscopic (SMM 0224.6; see §6.3| ) or photometric 
(SMM 1122.2) redshift, however the uncertainties are fairly large 
on the latter source. Two additional ambiguous/tentative sources 
(SMMs 0224.5 and 2152.2) have photometric redshifts for their 
most probable identification (i.e. the brightest MIPS emitter) in the 
SMG search radius. Therefore, for the remaining 16 SMGs (not 
including potential cluster members or those without any counter- 
part identifications), we must rely on more crude redshift estimates 
based on flux ratio trends. 

In principle, we can calculate a redshift for each radio identi- 
fication based on the strong evolution of the spectral index between 
850 \im and 1.4 GHz emission in AGN and star-forming galaxies. 
However, there is significant uncertainty in this technique (see fig. 
4 in Webb et al. 200 5 ), as it is highly sensitive to the dust temper- 
ature and to AGN contamination of the radio emission ( [Carilli &] 
[Yun|1999| ). We therefore do not use this redshift estimate for indi- 
vidual sources, but instead calculate a median redshift for the radio- 
identified sample (5 SMGs) of z > 2 (assuming a dust temperature 
of ^ 40 K), which is consistent with other submillimetre surveys 
with deep radio imaging (e.g. [Small et al.|2000{[Ivison et al.|2002{ 
[Chapman et al.|2003[[2005l ). 

We also investigate the redshift distribution of the SMG sam- 
ple through a comparison of 24 |j,m and 1.4 GHz fluxes ( [Ivison etaL] 
[2007| ) in the three cluster fields with radio and MIPS imaging (RCS 
0224, 1419, and 2319), as shown in Fig.[T4] The SMGs with pho- 
tometric (SMMs 0224.2, 0224.3, 0224.5; open circles) or spectro- 
scopic redshifts (SMM 0224.6; open circle with box) all are con- 
sistent with the track for Arp 220, supporting their identification 
and redshift. The two SMGs with robust mid-IR and radio detec- 
tions but no redshift estimate (SMMs 0224.4 and 2319.1; filled cir- 
cles) are also consistent with the Arp 220 track at their respective 
cluster redshift, and therefore may also be cluster members (but 
see §6.3[ for an alternative explanation regarding SMM 2319.1). 
The three SMGs of RCS 1419 with counterpart identifications but 
no radio detections (SMMs 1419.1, 1419.3, and 1419.4; dotted 
blue lines) are given a 3 a upper flux limit at 1.4 GHz and have 
S24fj.m/ SiAGHz ratios consistent with higher redshifts (z > 1) as- 
suming an Arp 220-like SED. Based on the lower limit of their flux 
ratio, we also plot these sources (filled blue circles with double ar- 
rows) interpolated onto the Arp 220 track. Although the remaining 
SMG (SMM 2319.2; dotted green line) lacks a radio detection and 
redshift estimate, it is not inconsistent with a lower redshift source 
based on the interpolated value (filled green circle with double ar- 
rows) to Arp 220 and could in fact be a cluster member. 

The submillimetre flux-density ratio provides another means 
of constraining redshifts (e.g. Hughes et al.[[1998] [Small et al.[ 
[T9991 [Fox et al.[[2002l [Coppin et al. 2008 ). In Fig. [15] we plot 
S 850 fim/S 450 fim as a function of redshift for Arp 220, M82, our 
entire SMG sample within the 450 \im field of view (22 SMGs to- 
tal), and four isothermal grey bodies with a range of dust tempera- 
tures (Td) and dust emissivities (e ocu^). Given the monotonic rise 
in the Arp 220 flux ratio with redshift, we plot the 3 a upper flux 
limit for 450 \im non-detections as a lower limit in redshift space 
(the ten cases of circles with arrows). SMGs with optical redshifts 
are shown as open circles (boxes represent spectroscopic redshifts), 
while those lacking redshift information (filled circles) are given in- 
terpolated values to the Arp 220 curve based on either flux ratios 
or limits. Ratios below the lowest value from Arp 220 are plotted 
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Figure 14. The ratio of S24nm/ Siaghz versus redshift for Arp 220 (solid 
line), Mrk231 (dotted Une; an AGN-dominated ULIRG), and for the ten 
SMGs with counterpart identifications in the three fields with radio and 
MIPS imaging. The four radio-identified SMGs (SMMs 0224.2, 0224.4, 
0224.6, and 2319.1) are plotted as open and filled circles without arrows; 
SMM 2319.3 is omitted since the 24 |j.m flux is severely blended. The re- 
maining six SMGs with MIPS emission but no radio detection are given a 
3 cr upper flux limit at 1.4 GHz and shown as circles with arrows. SMGs 
with (without) optical redshifts are plotted as open (filled) symbols; the 
spectroscopic redshift for SMM 0224.6 is highlighted by a box. The two 
larger open circles denote the only two SMGs with detections at both 24 \xm 
and 1.4 GHz and redshift estimates (SMMs 0224.2 and 0224.6). The four 
SMGs that are only detected at 24 |Ltm and lack redshift estimates (SMMs 
1419.1, 1419.3, 1419.4, and 2319.2) are plotted as horizontal dotted lines 
and at redshifts interpolated onto Arp 220 given their flux ratio limit; they 
are given both lower limits on their flux ratio and redshift (filled circles 
with double arrows). The dashed vertical lines represent the three cluster 
redshifts. Error bars on each point have been omitted to avoid crowding the 
plot; instead we show a representative error bar on the left side of the plot. 



arbitrarily as lower limits at z = 0.1. The plot clearly illustrates the 
sensitivity of redshift with the chosen SED for an individual SMG, 
and therefore, the colour ratio should only be used to estimate the 
redshift distribution of the entire sample (see also |Amblard et al.] 
|201 ). Assuming an Arp 220-like trend, our submillimetre flux ra- 
tio predicts that roughly over half our sample lies at z > 1 and the 
rest within z < 1. 

We note that for individual galaxies there is a fairly large scat- 
ter in the redshifts estimated from the techniques discussed above. 
Since the MIPS flux could suffer from confusion and the submil- 
limetre colour ratio is highly dependent on dust temperatures, the 
uncertainty on an individual redshift is typically ^0.5. However, 
we attempt to extract broad redshift trends for the sample; overall, 
we find that the combination of flux ratios between radio, MIPS, 
and 450 \im indicates that the majority of the our SMGs lie at 
z > 1. Even after accounting for individual redshift uncertainties, 
this implies that the bulk of our sources are in the background of 
the clusters. 



6.3 Galaxy-Galaxy Lensing 

Given the steep submillimetre source counts at the bright tail, re- 
cent studies have investigated the role of gravitational lensing in 
the SMG population (e.g .|Perrotta et al.'2003','Almaini et al.|2003| 
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Figure 15. The ratio of S'850/xm/*S'450/xm as a function for redshift for all 
SMGs within the 450 |LLm FOV. The soHd (dotted) line is the redshift track 
for Arp 220 (M82) and the dashed lines represent the range of redshifts for 
all seven clusters. We also plot isothermal grey bodies over a range of plau- 
sible dust temperatures and dust emissivities. The symbols are the same as 
in Fig. [14] Non-detections at 450 jitm are given 3 a upper flux limits (open 
and filled circles with arrows) and interpolated onto the Arp 220 track when 
possible. The two larger open circles denote the only two SMGs with de- 
tections at 450 M-m and redshift estimates (SMMs 0224.5 and 0224.6). Error 
bars on each point have been omitted to avoid crowding the plot; instead we 
show a representative error bar in the upper left corner. 



Austerma nn et al.||2009| [Negrello et al. 2010, Viei ra et al!] 
Aretxaga et al.|2 011). Theoretical models predict that ^ 30 



per cent of SMGs could be gravitationally lensed at fluxes above 



lOmJy at 850 ^tm fPerrotta et al.'2003'), with the fraction increas- 
ing at even brighter fluxes (Negrello et al. 2007). 

Current observational studies support this claim. For instance, 
[Almaini et al.| ( [2005^ ) found a positive angular correlation between 
SCUBA sources and low-redshift optical galaxies, with 20 — 25 
per cent of submillimetre sources associated with foreground over- 
densities within a 30 arcsec aperture. They argue that a possible 
cause of this correlation could be gravitational lensing. More re- 
cently, [Vieira et al.| ( [20 10] ) discovered 20 dusty sources above 
lOmJy at 1.4 mm within 87 deg^ with the South Pole Telescope 
(SPT), and postulated that they represent a population of gravi- 
tationally lensed, dust-dominated sources. Indeed, Negrello et aL] 
P010| ) claim that bright 500 |j,m-selected sources are almost all 
lensed. Moreover, [Austermann et al.]p00 9^ report a ~ 3 cr excess 
of SMGs compared to the blank-field population with the AzTEC 
camera in the Cosmological Evolution Survey (COSMOS), a 0.15 
deg^ field that includes a massive z = 0.73 cluster and an over- 
density of optical-IR galaxies at 2; < 1.1. They similarly calculate 
a spatial correlation between bright SMGs and low-redshift mass 
structure, attributable to lensing of the background submillimetre 
population by individual galaxies. [Aretxaga et al.| (2011) corrobo- 
rate this claim with a larger AzTEC survey of the COSMOS field 
(0.72 deg^) and find an excess of number counts compared to those 
of SHADES at 1 . 1 mm ("Auste rmann et al.|2010| ) due to moderate 
magnification levels from galaxy-galaxy and galaxy-group lensing. 

In light of these recent observations, we explore the possibility 
of galaxy-galaxy gravitational lensing within our sample between 
an optical cluster member and background SMG. Although lensing 
events are expected to be fairly common with the submillimetre 
population due to its high median redshift (z > 1) (Blain 1996), 
large-area surveys are essential to obtain a statistically significant 
sample of lensed SMGs since they have a low surface density at 
higher fluxes where they are expected to dominate over un-lensed 
SMGs ( NegreTT o et al.[[2010| ). Given the increased density of op- 
tical sources in cluster fields, there could be a larger correlation 
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in our sample. However, even a handful of gravitationally lensed 
SMGs in our relatively small sample, although difficult to statis- 
tically confirm, is noteworthy, as it is a possible complication for 
future observations of cluster fields. 

Again, we are biased towards our fields with radio and optical 
coverage, as accurate positions are necessary to differentiate the 
foreground optical lensing galaxy from the background SMG 
with radio emission. We find two possible cases of galaxy-galaxy 
lensing out of our five radio-identified SMGs. We note there are 
other cases of blended IRAC sources elongated in the direction 
of a possible background source (SMM 2152.1 and 2152.2), but 
they lack follow-up radio observations to verify the submillimetre 
position. Here we describe the counterpart identifications for the 
two possible anomalous lensed systems and include their postage 
stamps for visual reference (Fig. [16]) along with archival HST 
imaging (Fi^ 



17 1 to support this idea. 



SMM 0224.6. Counterpart emission is detected at mid-IR, 
radio, and optical wavelengths, all coinciding within ^ 2 arcsec of 
the 850 |j,m emission. Distant 450 \im emission at 3.1 a is detected 
11.4 arcsec away from the 850 \im source. The radio centroid lies 
between two bright optical sources and directly on top of diffuse 
HST emission (see Fig. [Tt]) and a faint i^-band detection (see 
[Webb et aL]|2005| ). The northern optical source has a foreground 
photometric redshift estimate of z = 0.351 =b 0.116, while the 
southern source has a spectroscopic redshift of z = 0.7832, 
consistent with the redshift of the cluster. 

Spectroscopic observations were obtained with the DEep 
Imaging Multi-Object Spectrograph (DEIMOS; Faber et al.|20"03] ) 
on the Keck II 10-m telescope as part of the Observations of Red- 
shift Evolution of Large Scale Environment survey (ORELSE; Lu- 
|bin et al.|2009| ). The diffuse HST emission, which is likely the sub- 
millimetre counterpart based on the radio position, was serendip- 
itously observed with DEIMOS and has a spectroscopic redshift 
of z = 1.0995. This system, therefore, represents a possible case 
of a background galaxy that may be lensed by either a foreground 
galaxy or cluster member. Assuming a singular isothermal sphere 
(SIS) profile and either a massive lens for the cluster member 
(z — 0.7832) or a lower mass (^ 60 per cent of the cluster member 
mass) for the z = 0.351 source results in magnifications of 1.3 
and ~ 1.7, respectively. This range of values, depending on the 
lens, is comparable to the possible magnification due to the cluster 
potential. As this is only one of two SMGs with a spectroscopic 
redshift in our sample, we briefly outline its properties. 

The spectrum of the submillimetre counterpart contains a 
faint continuum and a single emission line detection of strong 
[Oil] A3727A emission. The line flux of the [Oil] emission was 
measured using bandpass measurement techniques performed on 
the rest-frame spectrum, adopting the [Oil] bandpass of |Fisher| 
|et al.| (|1998|) an d a constant slit throughput of ujsut = 0.37 
( [Lemaux et al.||20T0] ). This measurement resulted in a line lu- 
minosity of L([OII]) = 2.75 X 10^^ ergs s~^, which translates 
to an [OII]-derived SFR of 3.85 Mq yr~^ using the relationship 
of |Kennicutt| (1998a). Combining this measur ement with the 
24^m-derived SFR of 123M0yr"^ from the [Chary & Elbaz 



^ Based on observations made with the NASA/ESA Hubble Space Tele- 
scope, and obtained from the Hubble Legacy Archive, which is a collabo- 
ration between the Space Telescope Science Institute (STScI/NASA), the 
Space Telescope European Coordinating Facility (ST-ECF/ESA) and the 
Canadian Astronomy Data Centre (CADC/NRC/CSA). 



( |2001| ) SED templates results in an E{B - V)s = 1.13. This 
value is consistent with other measurements of IR-luminous 
star-forming galaxies at similar redshifts (e.g. |Kawara et al.|2010| ), 
thus increasing our confidence that this galaxy is correctly matched 
to the submillimetre emission. As noted before in the literature 
(e.g. |Geach et al.||2009) , the optical star formation estimate is 
drastically lower than that from the mid-IR. 

SMM 2319.1. Radio, mid-IR, 450 \xm, and optical emission 
are all within 2.5 arcsec of the 850 |^m (5.5 a) source (see Fig. [16]), 
though the correct ID is ambiguous. The radio centroid lies slightly 
offset from the optical galaxy, as seen in the HST stamp (Fig.[T7]), 
and therefore might be a distinct source. Indeed, the IRAC 
emission looks like two blended sources: one coincident with the 
optical source, and the other associated with the radio source. The 
bright optical galaxy has a photometric redshift consistent with the 
cluster at z = 0.868 =b 0.139 . The IRAC colours further support 
membership as the source is consistent with the cluster redshift 
along the irregular starbursting galaxy track on the colour-colour 
diagram. The MIPS emission is closer to the optical source than 
radio centroid, but could be from both given the large PSF. A 
smoothed HST image (Fig. [TT]) reveals tentative evidence for 
faint optical emission at the radio position, which could indicate a 
background galaxy. This therefore might represent another case of 
galaxy-galaxy lensing between a cluster member in the outskirts 
and a background radio source. Conservatively assuming a redshift 
of ^ = 4, a fairly massive lens with a velocity dispersion of 
cr ~ 200 kms~^, and an SIS profile, we estimate a magnification 
of ^ 1.3. This modest magnification would be low compared to 
that from the cluster potential (^ 2.3) assuming the same source 
plane redshift of z — 4. We note, however, that the combination 
of the cluster mass distribution and galaxy-galaxy lensing can 
enhance the magnification further, as they do not add linearly. 



6.4 The Possible Origin of Excess Counts in Cluster Fields 

Contamination by both submillimetre cluster members and strong 
lensing by individual massive galaxies (either within the cluster or 
along the line-of-sight) might drive excess source counts in cluster 
fields, given the observed evolution of cluster galaxies and the re- 
cent detection of lensed SMGs, respectively. We describe each of 
these below. 

The growth and evolution of galaxy clusters, through the ac- 
cretion of galaxies, galaxy group infall, and cluster-cluster merg- 
ers, have drastic influences on populations within clusters and may 
temporarily spark extreme starbursts and fuel AGN activity ( Blainj 
let al.|1999a[|Owen et al.n 1999, Mi ller & Owen|2003| ). Moreover, 
since clusters accrete populations from the surrounding field, any 
evolution in cluster galaxies may trace trends in the field (e.g. |Tran| 
|inLp005| ). 

The first evidence for a rapid evolution of cluster galaxies 
emerged at optical wavelengths when [Butcher & Oemler] ( |1978| 
J 984 ) observed a larger fraction of blue star-forming galaxies in 
moderate redshift galaxy clusters {z ^ 0.5) compared to their lo- 
cal counterparts; this discovery has since been confirmed through 
many complementary tec hniques (e.g. [Dressier & Gunn|p^982[ 
[Couch & Sharples|[M7l [Ellingson et al.||2001| ). Local clusters, 
however, are dominated by passively evolving elliptical galaxies, 
suggesting that star-forming galaxies are transforming into the 
early-type population prevalent in clusters at z ^ ('PoggiantiJ 
|et al.jl999j ). Moreover, the field population, which feeds the infall, 
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Figure 16. 24 x 24 arcsec^ postage stamps for the two cases of possible galaxy-galaxy lensing: SMGs 0224.6 (top) and 2319.1 (bottom). The grey-scale 
images and symbols are the same as those in Fig. [73] 
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Figure 17. 10 x 10 arcsec^ HST postage stamps for the two cases of possible galaxy-galaxy lensing. The rightmost stamp is the HST stamp of SMM 2319.1 
smoothed with a 0.3 arcsec Gaussian to highlight seemingly faint optical emission at the radio centroid. The red contours show 1.4 GHz emission starting at 
45 |jJy and increasing in 10 |jJy steps. 



parallels this trend with an increasing star formation density with 
redshift (e.g. |Lilly et al.|1996[ [Ellis et al.|1996| ). 

Since the most extreme activity in field and cluster galaxies 
may be enshrouded by dust (e.g. [Due et alT]|2002| ), recent studies 
have transitioned to longer wavelengths in the infrared and submil- 
limetre to penetrate the heavy obscuration, and to more accurately 
measure the total star formation. Thus, the now well-established 
increase in the number of dust-obscured submillimetre luminous 
field galaxies with redshift ( [Smail et al.||1997| ) may similarly be 
mirrored in cluster environments. Indeed, evidence for enhanced 
star formation activity in moderate and higher redshift clusters has 
recently been reported at mid-IR wavelengths ( ,Fadda et al.|2000[ 



Geach et al.|2006[ |Bai et al.|2007[ [Saintonge et al.|2008[ [Hames] 

et al. 2009 ). If some of these systems reach ULIRG-level luminosi- 
ties they will be detectable with SCUBA at 850 |^m, and the trend 
should be even more apparent at the higher redshifts in our sample. 
Therefore, any submillimetre cluster member should be removed 
from the source counts for direct comparison to the blank-field. 

Given the recent detection of a population of gravitationally 
lensed dusty objects ( [Negrello et al.|2010[rVieira et al.|2010| ) and 
the correlation between observed SMGs and low-redshift galaxies 
( [Almaini et al. 12005) [Austermann et al.|2009| ), the excess of counts 
at higher fluxes could also be due to contamination from lensing 
effects. Although our sample only contains a few potential cases. 
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they could considerably affect the shape due to the steepness of the 
counts curve (Negrello et al.|2010| . 

To determine the potential origin of excess source counts 
in cluster fields, we have attempted to isolate these two possi- 
ble drivers. First, using the strongest member candidates (SMMs 
0224.2, 0224.3, and 2319.3), we investigate the effects of clus- 
ter member contamination in our source counts by removing these 
three galaxies from the calculation. The resulting cumulative counts 
(Fig. [18]) display minimal differences from the original counts. As 
members are only found within the lensing sample (but, again, this 
may be biased by the more extensive follow-up data), the source 
counts for the control cluster remain unaffected. 

We next exclude all possible cases of contamination from the 
source counts, including cluster members (SMMs 0224.2, 0224.3, 
and 2319.3) and possible examples of lensed SMGs (SMMs 0224.6 
and 2319.1). While this does bring down the cumulative counts 
of the lensing sample to blank-field levels, it is not sufficient to 
entirely remove the discrepancy between the two cluster samples. 
Further high-resolution observations could isolate lensed from un- 
lensed SMGs by providing clear positional information and possi- 
bly verify the source of discrepancy between the counts in clusters 
and blank fields, and moreover between the two cluster samples. 

We note that the other cluster surveys plotted ( [Cowie et al.| 
|2002[[Smail et al.|2002 ) have removed central cluster galaxies from 
their source counts analysis and corrected for strong lensing due to 
the cluster potential, yet still show an excess compared to the blank 
field. We do not attempt a lensing correction from the cluster po- 
tential, as the average lensing magnification is expected to be < 2 
beyond the central ^ 30 arcsec, even at the most extreme source 
plane redshifts. In our lensing sample, the two SMGs most centrally 
located in the clusters are SMMs 0224.6 and 2319.3, both of which 
have already been excluded in the right panel of Fig. [18] More- 
over, in the control group, only SMMs 1122.3 and 1326.1 (both 
within ^ 20 arcsec, projected) are likely to be significantly lensed 
(assuming source plane redshifts z > 2.5); removing these from 
the counts analysis would only amplify the discrepancy between 
the two samples. Future cluster surveys should ensure all possible 
members are accounted for and possible lensing effects due to in- 
dividual galaxies and the cluster potential are noted. 

The fact that all of our likely cluster members originate from 
lensing clusters might be indicative of dynamically unrelaxed states 
or physical mechanisms (e.g. mergers of cluster galaxies) that en- 
hance activity and induce star formation, while temporarily boost- 
ing the lensing cross section. Indeed, RCS 0224 represents one 
of highest known redshift systems with such bright and abundant 
lensed sources ( |Gladders et al.|2002j |2003| , and contains at least 
three optical arcs, the largest of which originates from a back- 
ground source at 2; = 4.88. Moreover, preliminary measurements 
of the dynamical state from 117 member galaxies observed in the 
ORELSE survey show evidence for an unrelaxed cluster with sig- 
nificant substructure. RCS 2319, a rare and distant supercluster, 
comprises three spectroscopically confirmed bound clusters, sep- 
arated by a projected distance of < 3Mpc, which are expected 
to merge by z ^ 0.5 to form a massive structure with mass 
^ 10^^ Mq jGilbank et al.|2008|. The cluster wit h SCUBA cov- 
erage has three optical arcs ( [Gladder s et al.||200'3] ), one of which 
belongs to a high-redshift Lyman break galaxy at z = 3.86. 

However, the lensing clusters seem to be fairly relaxed sys- 
tems in terms of their X-ray properties; in particular, RCS 1419 is 
compact with very regular X-ray morphology. On the other hand, 
RCS 2318 of the control sample shows significant substructure at 
X-ray wavelengths (Hicks et al.|2008j . Perhaps then, the dynamical 



distinction between the samples is that of a relaxed cluster with a 
sub- structured collapsed core (i.e. the super-lenses) versus a wholly 
unrelaxed cluster. The former may therefore be a signpost for both 
enhanced lensing due to a boosted central lensing cross-section and 
increased star formation in the form of excess SMGs, for example. 



7 CONCLUSIONS 

We have presented a study of seven galaxy clusters between 
0.640 < z < 1.045, three of which display strong gravitational 
arcs. Within these fields, which cover a total of 47 arcmin^, we 
have detected 26 SMGs and SMG candidates at 850 M,m with S/N 
^3.0. We summarize our results as follows: 

(i) As with other cluster surveys, we see an excess of SMG 
source counts compared to the blank field; however, this trend is 
limited to our strong lensing sample. The control sample actually 
shows a slight deficit of counts, though consistent with the re- 
sults from the LABOCA ECDFS Submillimetre Survey (LESS); 
the ECDFS has been shown to be an underabundant field at many 
wavelengths (WeiB et al. 2009 and references therein). There is a 
2.7 a excess of 850 |j,m source counts in the lensing sample com- 
pared to the control clusters. 

(ii) We utilize mid-IR, radio, optical, and X-ray catalogues to 
determine counterpart emission to the SMGs, thereby improving 
the SMG astrometry and, in some cases, providing redshift esti- 
mates. We classify our identifications as secure, ambiguous, or ten- 
tative, and identify 12, 6, and 5 SMGs for each classification, re- 
spectively. Three SMGs lack any plausible counterpart emission 
and are therefore indicative of either spurious sources or high- 
redshift galaxies. 

(iii) Employing infrared colour criteria and locating optical 
counterparts with photometric redshifts allows us to identify three 
likely cluster members: SMMs 0224.2, 0224.3, and 2319.3. We cal- 
culate infrared luminosities of ~ 10^^ Lq, consistent with ULIRG 
levels. Assuming no AGN contamination, all three members have 
extreme star formation rates of ^ 700-1000 Mq yr~^, but we ar- 
gue that the most luminous of the three probably contains an AGN. 
The extreme luminosities and star formation rates of these ULIRG- 
level sources are not typical in cluster environments, which are 
usually populated with gas-poor spirals or old elliptical galaxies 
( Sand ers & MirabeT| ri996 ) at low redshift. Therefore, these cluster 
members may reflect the increased infall and activity that is seen in 
clusters at higher redshift with the Butcher- Oemler effect and sim- 
ilar mid-IR studies. All three likely cluster members derive from 
strong lensing clusters that display evidence of additional cluster- 
wide activity: the RCS 2319 supercluster comprises three merging 
systems that give rise to disturbed galaxy morphologies; it is con- 
ceivable that RCS 0224 might be undergoing similar activity, al- 
though thus far, current measurements in the ORELSE survey only 
hint at enhanced substructure. 

(iv) We observe two possible cases (SMMs 0224.6 and 2319.1) 
of galaxy-galaxy lensing between an optical foreground galaxy 
(or cluster member) and a background SMG. Follow-up high- 
resolution observations or the detection of CO emission lines are 
essential to verify this claim, but beyond the scope of our current 
work. 

(v) While there is some evidence for cluster member contami- 
nation in the source counts, it is not sufficient to remove the excess 
from the lensing fields. We note, however, that the effect of cluster 
contamination is only a lower limit as we are biased by the quantity 
of the multi- wavelength follow-up. Galaxy-galaxy lensing can also 



© 2011 RAS, MNRAS 000,[l]|28] 



24 A.G. Noble et al 



10^ 



10" : 



10^ 



10^ 



10^ 



10° 



o This work - Lensing Clusters 
O This work - Control Clusters 

Small et al. 2002 

Cowie et al. 2002 
)K Webb et al. 2003 - CUDSS 

Scott et al. 2002 - 8 mJy Survey 
X Coppin et al. 2006 - SHADES 
o Weiss et al. 2009 - LESS 



4' 



<> 



0.1 



1.0 

Flux (mJy) 



10.0 



10^ 



10" 



10^ 



10^ 



10^ 



10« 



<> 



o This work - Lensing Clusters 
O This work - Control Clusters 

Small et al. 2002 

Cowie et al. 2002 
)K Webb et al. 2003 - CUDSS 

Scott et al. 2002 - 8 mJy Survey 
X Coppin et al. 2006 - SHADES 
o Weiss et al. 2009 - LESS 



0.1 



1.0 

Flux (mJy) 



10.0 



Figure 18. [Left] Cumulative 850 |Lim source counts, excluding the three likely cluster members. Our results are denoted by the larger symbols, where the 
red circles represent counts from the strong lensing clusters and the blue diamonds represent counts for the control group. [Right] Cumulative 850 |^m source 
counts, excluding the three cluster members and two possible cases of galaxy-galaxy lensing. 



contribute to source count contamination at higher fluxes, but the 
effect on our counts must be verified through follow-up observa- 
tions on the possible lensed systems. We suggest that the dynamical 
state of the cluster might account for any additional source count 
discrepancies between the two cluster samples, as a sub- structured 
core, of an otherwise relaxed cluster, could boost central lensing 
efficiency while triggering extreme starbursts and driving the for- 
mation of ULIRGs within cluster fields. 

(vi) While we have not unambiguously accounted for the ex- 
cess source counts between the two cluster samples, we have high- 
lighted possible problems for cluster surveys, including increased 
cluster member contamination at higher redshifts and potentially 
lensed SlVlGs. Further work is required to fully isolate the driver of 
contamination, such as high-resolution Submillimetre Array imag- 
ing, CO emission line observations, and homogeneous radio cover- 
age. Although this study suffers from both small number statistics 
and difficulties associated with the SIMG counterpart identification 
process, it still has important ramifications for intrinsic submillime- 
tre source counts and SZ cosmological studies that are concerned 
with point source contamination. We therefore urge caution for fu- 
ture large-scale cluster surveys, for instance with the Herschel tele- 
scope and the forthcoming SCUBA-2 instrument, and stress the im- 
portance of accurate lens models and proper subtraction of cluster 
members for counts analyses. 
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APPENDIX A: A DISCUSSION OF INDIVIDUAL 
SOURCES 

Here we provide a discussion of each submillimetre source detected 
in the survey. We reiterate that individual sources detected below 
3.5 cr have a high probability of being spurious, making counter- 
part identification rather uncertain. Nevertheless, these sources are 
included in the discussion for completeness. Table [5] summarizes 
our multi- wavelength counterpart analysis and Figures [A 1 1 A5 | pro- 
vide visual support for each description. 



Al Secure counterparts 



6.1.1 



SMM 0224.2. See J 
SMM 0224.3. See | 
SMM 0224.4. We find secure MIPS, radio, IRAC and 450 ^im 



6.1.1 



counterparts within the lOarcsec search radius. The 450 ^im emis- 
sion (3.3 cr) is isolated, offset in declination from the other three 
sources. The MIPS source is 5.6arcsec away from the 850|j,m 
emission but within 0.5 arcsec of the radio and IRAC positions. 
Based on positional offsets to the MIPS emission and infrared 
colours, the more northern of the two blended IRAC sources 
(shown in the 3.6 ^tm grayscale image) is the more likely SMG 
counterpart. There is also a faint i^^-band detection at the radio cen- 
troid with extremely red colours (see |Webb et al.| ( [2005] )), consistent 
with a high-redshift galaxy. It is too faint to for sufficient detection 
in optical bands, and we must therefore rely on the redshift deduced 
from the radio flux, z = 2.56, making the SMG consistent with a 
background galaxy. HST imaging (see Fig. |A5| l reveals two optical 
sources a few arcseconds south of the radio centroid (also seen in 
the z' image). The western source is likely a star, while the eastern 
detection could be a foreground galaxy acting as a lens to the back- 
ground radio source, but further observations, beyond the scope of 
this paper, are required. 
SMM 0224.6. See ^3] 

SMM 1122.2. Within 10 arcsec, only one 24 |j,m source and one 
optical source are detected. The MIPS source is offset by 0.7 arcsec 
and therefore provides a secure identification. The optical source 
lacks a reliable photometric redshift (z = 1.073 =b 0.803), so we 
cannot determine cluster membership for this SMG, but HST imag- 
ing (see Fig. |A5| l reveals a faint compact source. Of the four IRAC 
detections, only one is significant at 3.6 |j,m and coincides with the 
MIPS and optical emission. This source lies within the likely SMG 
counterpart region of the colour-magnitude diagram (z > 1.5) and 
within the colour-colour boundary for SMGs. 
SMM 1122.3. An extremely bright MIPS source (1096.4 |aJy) lies 
7.5 arcsec away with coincident optical and IRAC sources. Slightly 
further out (9.9 arcsec from the SMG) is 450 \im emission detected 
at a S/N of 5.63. The optical counterpart lacks a photometric esti- 
mate, but the IRAC colours suggest this SMGs lies at z < 1.5. The 
HST image (Fig. |A5| l shows a likely foreground galaxy, possibly a 
face-on spiral with a central nucleus. We therefore conclude this is 
a secure ID, and most likely at lower redshift. 
SMM 1326.1. MIPS, IRAC, and 450 |^m (6.2 a) emission all ex- 
ist within 5 arcsec of the submillimetre source, offering a robust 
counterpart identification, but there is no visible optical compan- 
ion. The IRAC source falls within the colour-colour boundary from 
|Yun et al^ ( |2008| ) and within the colour-magnitude selection crite- 
ria for SMG counterparts proposed by |Hainline et al.| (|2009 ), and 
therefore likely at z > 1.5, which would also explain the lack of an 
optical detection. 

SMM 1326.2. The total 24 |j,m flux density within this search radius 
is significantly dominated by one MIPS source (249 pJy). We also 
detect X-ray emission at 5.6 a, coincident (offset by 1.1 arcsec) 
with the brighter MIPS source. As we do not have any photometric 
redshifts for this entire cluster field, we cannot determine cluster 
membership, but there is a faint optical counterpart coincident with 
the secure ID. The two IRAC sources associated with each MIPS 
object are strongly detected in every channel except at 8.0 |^m, so 
we can only place an upper limit on the Sg.o fim/ S4.5 fim colour. 
However, the redder S^.s ^m/ S^.q colour of the IRAC source 
associated with the X-ray and brighter MIPS emission places the 
object in the likely SMG counterpart region, making this a secure 
identification, likely at high redshift. 

SMM 1326.3. We find a single MIPS source (coincident with IRAC 
and optical emission) with a 2.6 arcsec offset from the SMG, pro- 
viding improved positional accuracy for the submillimetre emis- 
sion. The associated IRAC source secures the identification as it 
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falls within the colour-magnitude selection criteria given by |Hain-| 
nine et al.l2009l 

SMM 1419.1. A single faint MIPS source (63 |aJy) is located 
4.8 arc sec away and offers better astrometry for the submillime- 
tre source. The IRAC source associated with the MIPS emission 
is weakly detected, so we can only estimate a limit for its infrared 
colours. While the other IRAC source in the field is brighter in all 
channels (except at 8.0 |^m ), its infrared colours are inconsistent 
with an SMG, so we classify the MIPS/IRAC detection as the more 
likely counterpart. A smoothed HST image (see Fig. |A5| ) reveals a 
two possible optical counterparts: one just southwest of the MIPS 
centroid, and another fainter source aligned with the IRAC posi- 
tion. We note that there are no photometric redshifts available for 
this entire cluster. 

SMM 1419.4. We detect a 111 |j,Jy MIPS source within the search 
area, altering the submillimetre astrometry by 6.4 arcsec. There is 
also a coincident IRAC source, but weakly detected at 5.8 \im and 
without a detection at 8 ^tm, preventing a reliable colour analysis. 
The HST image (Fig. |A5| l shows a slight extension in the optical 
emission. As there is only one MIPS source in the entire search 
radius, we classify this as a secure identification. 
SMM 2319.2. We locate 24 |^m emission at 2.5 arcsec from the 
SMG. Corresponding 450 p,m emission (4.4 a) lies 1.1 arcsec be- 
yond the MIPS source, yielding a secure detection for this 7.5 a 
SMG. The associated IRAC source to the MIPS emission falls 
within the SMG counterpart boundaries on the colour-magnitude 
diagram (which typically indicates z > 1.5) and colour-colour dia- 
gram. Besides the star to the northwest, there is only a faint optical 
source (see Fig. |A5| l slightly offset from the IRAC position, possi- 
bly indicative of a high-redshift source. 

A2 Ambiguous counterparts 

SMM 1122.1. We discover two MIPS sources, each with a coin- 
cident optical and IRAC source, and both within 5 arcsec of the 
850 i^m emission, but this identification is complicated by a third 
MIPS/IRAC detection, blended with the northern MIPS source. 
Based on IRAC colours (although somewhat unreliable due to the 
confusion), the central IRAC source is the most likely counterpart. 
Deep HST imaging reveals faint optical emission between the two 
bright optical sources; this source becomes even more obvious af- 
ter smoothing the image with a 0.3 arcsec Gaussian (see Fig. |A5| l, 
and is probably a high-redshift source. Photometric redshifts for 
the two bright optical sources place them in the foreground at 
z = 0.282 ± 0.023 (northern) and z = 0.30 ± 0.032 (southern). 
We therefore classify this ID as ambiguous, but a possible exam- 
ple of galaxy-galaxy lensing between the northern foreground op- 
tical source and the central IRAC/MIPS/optical emission (our most 
likely SMG ID). 

SMM 1326.4. Confused MIPS sources populate the area surround- 
ing this SMG, complicating the ID. However, there is only one dis- 
tinct IRAC source with coincident X-ray emission as well. The red 
S5. 8 fim/ S3.6 ij.m colour placcs this source in a SMG region of the 
colour-magnitude plot (z > 1.5), supporting the MIPS and X-ray 
emission as the most likely identification. The infrared colours are 
within the colour-colour boundary, specifically along the power- 
law AGN region. Furthermore, the MIPS-to-IRAC colours are in- 
dicative of an AGN component with S2A^jim/ Ss ^im — 1-7 and 
Sg.o fim/ fim = 3.1 (Iv ison et al.||2004j ). Due to the blended 
MIPS emission, we classify this ID as ambiguous, but most likely 
a background source as it lacks optical emission. 
SMM 2152.1. Although we find two MIPS sources within the 



10 arcsec search radius, one is clearly a few sources blended to- 
gether. While slightly less confused, the IRAC emission can also 
not be separated into distinct sources for reliable colour estimates. 
The sole IRAC source associated with the fainter MIPS emission 
in the northeast has colours consistent with an elliptical galaxy, 
making it an unlikely SMG candidate. HST imaging (see Fig. \A5\ 
reveals faint emission to the southwest of the dominating fore- 
ground object. This therefore might be another case of galaxy- 
galaxy lensing, making the most likely submillimetre ID the back- 
ground lensed galaxy, but still rather ambiguous. 
SMM 2152.2. Severely blended MIPS emission is coincident with 
optical and blended IRAC sources. The bright optical source is a 
foreground galaxy (z = 0.301 =b 0.030), but a fainter source evi- 
dent in the HST image (Fig. |A5| l has a more uncertain redshift of 
0.827 =1= 0.468, which tentatively agrees with the cluster redshift. 
The IRAC and MIPS emission are too confused for reliable pho- 
tometry. We classify this as a possible galaxy-galaxy lensing pair, 
with the background source as the more likely SMG counterpart, 
but the ID remains ambiguous. 
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A3 Tentative counterparts 

SMM 0224.5. Three MIPS sources are detected within the 
10 arcsec search radius of this SMG. The brightest source (225 |jJy) 
dominates the other two (both are < 90 /xJy) by a factor ^2.5 and 
we therefore adopt this source as the primary contributer to the 
850 i^m emission. It is offset in right ascension from the SMG po- 
sition by 7.7 arcsec and in proximity of two optical sources. The 
closest (0.8 arcsec away) has a photometric redshift estimate of 
z = 0.943 =b 0.074, likely indicating a background galaxy, but 
consistent with the cluster redshift at just over the 2 a level. The 
corresponding IRAC source lies within error of SMG counterpart 
region on the colour-colour plot. Emission at 450 |^m also exists, 
closest to a weaker 24 |^m source, but we emphasize the high po- 
sitional uncertainty of such a low 450 \xm detection (2.8 a) and as- 
sume it too is associated with the brightest MIPS source. Due to the 
multiple MIPS detections, we label this identification as tentative, 
but place a higher counterpart likelihood on the brightest 24 \im 
source (most eastern). Given its low S/N at 850 |j,m (3.4 a), it could 
also be a blend of multiple submillimetre sources. 
SMM 1419.3. We find two MIPS sources, each with an IRAC de- 
tection, and 450 ^tm emission (at 2.7 a) neighbouringing the SMG. 
The MIPS sources are of roughly equal flux densities (^ 55 /xJy), 
but one source is closer to both the 450 and 850 \xm emission. 
The IRAC source associated with the closer MIPS position also 
has infrared colours more appropriate for an SMG. The source is 
firmly within the colour-magnitude boundary and situated along 
the power-law AGN track in the colour-colour diagram. The HST 
image (Fig. |A5| l reveals a compact optical object associated with 
the MIPS/IRAC emission, making this a likely, but tentative coun- 
terpart. 

SMM 2318.1. We detect strong 450 \im emission (7.1 a) coinci- 
dent (2.2 arcsec offset) with this extremely bright 850 ^im SMG 
(14.9 mJy with a S/N of 17.0). There is also a tentative optical 
counterpart close to the 450 i^m emission in the HST image with a 
disturbed morphology (see Fig. |A5| ). Given the steep drop in SMG 
source counts at flux levels greater than ^ lOmJy, such a bright 
source is fairly unique and could represent a strongly gravitational 
lensed SMG. Alternatively, it might be multiple detections of dis- 
tinct objects; indeed, the 450 |j,m emission is fairly elongated (see 
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Fig. [2]) and possibly double-peaked. Further high-resolution obser- 
vations are needed to verify either possibility. 
SMM 2318.2. Emission at 450 \im (3.0 a) lies within 2.2 arcsec of 
the submillimetre 850 ^tm source (5.6 a). The HST image (Fig. |A5| ) 
displays a few possible counterparts. 

SMM 2318.3. This SMG has 450 ^im emission, but is detected just 
within the search radius for 450 ^tm sources at 12.0 arcsec. HST 
imaging (see Fig. |A5| l reveals several tentative optical counterparts, 
but this SMG lacks any secure identification. 

A4 Insecure counterparts or possible spurious sources 

SMM 0224.1. The robust 850 \im detection of this source (6.1 a) 
means the assigned search radius is fairly small (7.4 arcsec). This 
source only has candidate counterpart emission at 450|j,m (3.7 a 
at 4.1 arcsec away) and two possible optical counterparts. Due to 
the high probability of detecting optical sources within this search 
area, we are unable to assign it a unique identification without a 
supporting radio or MIPS detection. We note, however, that there 
are radio and mid-IR sources (each with coincident IRAC and op- 
tical emission) just beyond the search radius with offsets of 9.9 and 
8.0 arcsec, respectively, but on opposite sides of the SMG. Given 
that the search radius is a 95 per cent confidence interval, it is pos- 
sible, though unlikely, that one of these objects is the counterpart. 
For now, however, we conclude that this source has no secure coun- 
terpart, perhaps due to a high redshift. 

SMM 1419.2. Weakly detected at 850 \im (3.0 a), this source only 
has counterpart emission at IRAC and optical wavelengths (but 
no photometric redshifts available). The brightest IRAC source at 
3.6 \im is an unlikely SMG candidate based on its infrared colours, 
and most likely a star based on the HST image (Fig. |A5| ). The best 
suited IRAC source is the most southern detection and is located 
within both appropriate SMG boundaries on the colour-colour and 
colour-magnitude diagrams. There is an HST source nearby that 
might be associated with the faint IRAC detection, however, with- 
out an associated MIPS source, we classify this identification as 
insecure. 

SMM 1419.5. This source evades counterpart emission at all other 
wavelengths, making it a possible candidate for one of the spurious 
sources. We note there is MIPS, IRAC and faint radio emission 
in the southwest corner, but beyond the assigned search radius of 
8.0 arcsec. Given the high detection at 850|j,m (6 a), this source 
could instead be a high redshift SMG, and therefore not sensitive 
to counterpart emission. 
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Figure Al. 24 x 24 arcsec^ postage stamps for SMGs with secure counterparts. The grey-scale images are optical z' band, IRAC 3.6 |xm, and MIPS 24 |xm, 
from left to right, where north points up and east is to the right. Each image is centred on the SMG position, denoted with a large blue cross. The size of 
the cross represents the search radius for MIPS counterparts. The rest of the symbols are as follows: an orange cross is 450 |j.m emission, a yellow square is 
a MIPS source, magenta triangles signify IRAC detections, the red contours show radio emission, purple diamonds show optical sources with photometric 
redshifts, and a green 'x' denotes an X-ray detection. Note that the radio contours cover the entire field, but the symbols are only shown for sources within the 
search radius. The radio contours start at 45 |Ljy and increase in 10 |jJy steps up to 105 |jJy, after which they increase by 200 [jJy. 
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